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ABSTRACT 


The use of remote sensing data for surface mapping is very common and is being 
carried out on routine basis. The IRS-IC LISS-III false colour composite image of 
Kanpur area shows a large channel migration pattern through neck cut off along Ganga 
river. This cut off spreads over an area of about 150 km^ showing that the Ganga river 
has migrated 15 km in the east from the present position during past several years. In 
order to understand the migration pattern, detailed image processing of data have been 
carried out to delineate the fractures and faults in the cut off area. Efforts have been made 
to use different edge operators: Sobel, Laplacian, and High pass filters (HPF). It has been 
found that only Sobel and High pass filters give good results. Otsu thresholding technique 
has been applied before and after merging, and binary images have been generated, 
which contain two pixel value ‘0’ and ‘255’, higher pixels value ‘255’ represent an edge. 
The thresholding technique has provided the ability for knowing the fault location. To 
improve the spatial resolution as well as to improve the visual interpretation, and to see 
the effect of merging on thresholding technique, various merging methods: HPF, IHS 
(Hue Intensity Saturation), P+Xs, and Price algorithm have been applied. The HPF 
method found to give better results due to different time for PAN and multi-spectral data. 
The HPF and Price methods have been used successfully to remove cloud cover on high 
resolution images. I’hc normalized difference vegetation index (NDVI) image of the area 
has been generated which shows significant difference across the expected fault. The 
analysis of soil across and quality of water samples in the area further support the 
presence of subsurface fault in the area. We have carried out 14 resistivity soundings in 
the area; the sounding results show drastic change across a profile crossing the large cut 
off. The integration of remote sensing data, surface and subsurface information deduced 
from resistivity soundings has been discussed in the Ganga river migration near Kanpur. 
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Chapter I 


INTRODUCTION 


1.1 GENERAL 

Remote sensing has proved as an important tool to provide basic information for 
the evaluation of natural resources. Remote sensing technology has played a major role in 
mapping of the Earth’s resources using high-resolution data. Presently, digital analysis of 
remote sensing data has become less expensive and more effective in mapping resources 
for various purposes. Conventional geological and geophysical surveys require much 
time and huge expenditure for mapping of inaccessible mountainous regions. Missing of 
rock outcrops and limited view of geological features, as well as manual methods of data 
collection and processing are major difficulties in the geological information system 
(GIS). Using remote sensing data a continuous record of information about ground 
features is obtained. High resolution data and synoptic nature of the information are 
useful in the interpretation of geological data. However, the degree of geological 
information processed from remote sensing data remains far behind that processed 
through geological and geophysical surveys. This is due to the limitation of present 
sensors in identifying rock formation and surface features. However, efforts are going on 
to improve methodologies for the interpretation of remote sensing data. Ground 
geological and geophysical surveys only confirm the geological features. The integration 
of geological, geophysical and remote sensing system will help to interpret the geological 
information. Various geological themes for which, both airborne as well as satellite borne 
data have been used to map the lineaments, lithological, structural, and 
geomorphological, mineral resources, and ground water targeting, engineering geological 
applications, river basin studies, plate tectonics and earthquakes studies, and 
environmental assessment. 

Remote sensing data has been used widely for different purposes in the 
applications of sciences and engineering. Image processing involves the manipulation 
and interpretation of digital images with the aid of a computer. It has huge numbers of 
applications, such as, image transmission and storage for business applications, medical 
processing, radar, industrial applications, sonar and acoustic image processing. Images 
acquired by different kinds of satellites with different bands are useful in tracking of earth 
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resources; geological mapping, such as lithological mapping, earthquake and volcanic 
monitoring; agricultural crops, urban growth, and weather; flood and fibre control; and 
many other environmental applications. 

Recently, there has been a remarkable increase in the number of available satellite 
and airborne sensors. Each sensor has its mission and applications. For many applications 
the combined information from multiple sensors provides more comprehensive 
information by collecting a wide diversity of sensed wavelengths, spatial resolutions, and 
look angles. This led to develop a new technique called image merging. 

The merging of multi-sensor and multi-resolution data is becoming a widely used 
procedure because of the complementary nature (different spatial and spectral resolution) 
of the various data sets. After image merging, we have high spectral resolution data set 
with high spatial resolution. Several merging techniques have been developed in the past 
to merge multi-sensor and multi-resolution data. Ideally, the method for merging data sets 
with high spatial and spectral resolution should not distort the spectral characteristics of 
the high spectral resolution data. Distorting the spectral characteristics of the data is not 
good for calibration purposes between original merged images. The main aim of data 
merging is to improve the information content of the images for visual interpretation and 
also improve the spatial resolution. 

1.2 OBJECTIVES OF THE PRESENT THESIS: 

i) To study the effect of thresholding technique on edge extraction. 

ii) To map the surface features using remote sensing data. 

iii) To map the surface and subsurface structure using remote sensing data. 

iv) To study the resistivity variations in the area. 

v) I'o find a remarkable explanation for the migration pattern of Ganga river in 
Bithoor area. 

The broad objective of the present research is to know the subsurface structure by 
using remote sensing data, and to find the effect of merging and thresholding techniques 
on edge extraction. Finally, to integrate remote sensing data with the resistivity variations 
obtained from conventional resistivity method for the same area. 



1.3 OVERVIEW OF STUDY AREA AND DATA RESOURCES 

The work presented in the present thesis has been conducted in agricultural area 
of Bithoor-Unnao Kanpur, Uttar Pradesh. The study area of Bithoor is located between 
North latitudes 26° 30' 30" - 26° 45' 30" and East longitudes 80° 15' 30" - 80° 30' 30". The 
Kanpur region has a flat topography with large part of agricultural area. Ganga river 
shows migration pattern through cut-off, which has paramount importance in 
understanding the tectonics, and subsurface dynamics, and migration pattern of Ganga 
river. The prominent meander cut-off spread over an area of about 150 km^ along the 
Ganga river in the Kanpur region. The subsurface morphology of the basin and tectonics 
must have played an important role in the sedimentation pattern of the region. 

Details of the satellite data products used in the present study are given in the 
Table '1.1. These data have been acquired from National Remote Sensing Agency 
(NRSA), Hyderabad. For reference, map products have also been used. Topographic 
maps have been obtained from Survey of India (SOI). Various images and maps used for 
the study arc listed below: 

1. Topographic map no. 63 B/ 6 on a scale 1: 50,000 to know various 
details of the area. 

2. IRS-1 C LISS-III multi-spectral satellite data acquired on February 
26^ 1997. 

3. IRS- 1C LISS-III panchromatic satellite data acquired on December 
16^ 1996. 

4. OYO Resistivity meter for carrying resistivity fieldwork. 

5. GPS has been used in the field to determined longitudinal and 
latitudinal coordinates. 

Table 1.1 Satellite data products used for the present study 
(from IRS- 1C user manual, 1995) 

Satellite f Band [ Resolution Size Wavelength Spectral region Path-Row 

(m) (pixels) (lim) 

"iRS'-Tc ' BarKlT 23J5 1057x987 0.52-0.59 Green 

LISS-in “Band! 2l75 1057x987 0.62-0.68 R^d 99 - 53 

Band a | 23.75 1057x987 | 0.77-0.86 Near infrared 

0.50-0.75 








1.4 SOFTWARE DETAILS 

Programming m C language has been developed for different merging methods 
and for thresholdmg techniques and for different edge operators to extract edges from 
images. Standard package ‘TDRISI for windows” which is an mtegrated package for 
image processing and GIS have also been used. This package was used for the 
registration of images, and for on screen digitization of the reference map. Resixp 
package is used for the mteipretation of resistivity data, which was collected through 
fieldwork. 

1.5 ORGANIZATION OF THE THESIS 

The thesis has been organized in five chapters. The present chapter provides a 
brief discussion on need of remote sensing data (data merging and thresholding 
technique) and geophysical data for determining the subsurface information, followed by 
details of study area and data resources. Chapter II deals with the geological and 
geophysical information of the study area. Chapter IE summarizes the earlier research 
work and theory of image processing and electrical resistivity method. This Chapter also 
presents in detail the theoretical backgroxind of data merging and edge extraction with the 
emphasis to various techniques implemented in the present thesis. Chapter IV presents 
various experiments and results and shows the correlation of both remote sensing and 
geophysical data. Finally, Chapter V presents conclusions of the present study. 



Chapter n 


Geological and Geophysical Information 
About the Study Area 


2.1 GENERAL 

Extensive geological and geophysical surveys have been carried out in Ganga 
basin since the last century. The Oil and Natural Gas Commission particularly intensified 
the investigations during the last three decades for hydrocarbon exploration. These 
surveys have brought out a clear picture of the tectonic framework. Oldham (1883) and 
Meddlemiss (1890) conducted the systematic geological surveys in the Himalayan 
foothills in the last centuiy. Wadia (1928), Auden (1934), and several other geologists of 
the Geological Survey of India have mapped the Tertiary sediments of foothills region. 
With the formation of Oil and Natural Gas Commission in 1956, efforts were intensified 
to carry out detailed geological mapping of the entire Tertiary belt. The important 
compilations have been made by Sahni and Mathur (1964), Karunakaram and Rao 
(1976), Rao (1974), Raiverman et aL (1976) highlighting the contribution of the ONGC. 
This area has also been covered by geophysical surveys (Das et al 1976), Balakrishnan 
and Chaudhury, (1976). The structure and stratigraphy of the alluvium covered Indo- 
Gangatic plain were subjects of speculation for a long period by many authors like Suess 
(1904), Burrard (1915), Oldham (1917), Wadia (1975), and Krishnan (1964). During the 
last three decades the geophysical surveys were carried out by the ONGC have brought 
out a clear picture of basement configuration of the area. Agrees (1957) has carried out 
aeromagnetic surveys over the Ganga valley. Systematic geophysical surveys have been 
carried out by ONGC (Chudhury, 1975, Ramachandra Rao, 1973, Balakrishnan and 
Chudhury, 1976, Aditya et aL, 1976). This chapter briefly discusses the geological and 
geophysical background of the Ganga basin. 

2.2 GEOLOGY OF THE GANGA PLAIN 

2.2.1 General Information of Ganga Basin 

Ganga plain of Uttar Pradesh, forming a part of the largest alluvial plain in the 
world, the Indus-Ganga-Brahmaputra plain, lies between the Peninsular Shield in the 



south and Himalaya in the north. It covers an area of about 2.23 laJch sq. Km. in U.P. 
Many theories have been put forward regarding the origin of this vast depression often 
referred to as “Indo-Gangetic Trough”. (Suess, 1893, Burrard, 1915, 1932, Oldham, 
1917, 1923, Wadia, 1938, 1966). However, it is now widely accepted that the Foreland 
basin (Le Fort, 1975) comprising Ganga plain, was formed due to continent-continent 
collision of Indian plate with Asian plate resulting in underthrustmg of the Indian plate, a 
phenomenon still in progress. The Ganga plain in U.P. lies between Yamuna in the west 
and 60 m topographic contour in the east. In the north, the Himalaya demarcates it. The 
Peninsular Shield marks the southern boundary. It comprises Upper G anga plain, 
Bundelkhand and Vmdhyachal-Baghelkhand regions (Singh, 1971). It is the oldest 
Quaternary platform formed by the paleo-drainage system; relicts and trails are still seen 
over it as left over features. The drainage network of the area comprises three major snow 
fed Himalayan rivers namely Ganga, Yamuna and Ghaghara The Ganga river being the 
principal trunk channel, is the recipient of all water lines in the region. The most 
conspicuous feature of the drainage configuration is the regional parallelism between 
courses of the rivers. This seemingly featureless plain lacks topographic prominence and 
the river bluffs, levees, and dead arms of the river channels and the river channel break 
the monotony of the physical landscape. The Trans-Yamuna area exhibits a distinct 
topography represented by deep valleys separated by sharp spins and buttresses due to 
ravine landscape topography along rivers and streams. This feature is not very 
uncommon with other streams of the rest of the area but hardly comparable in magnitude 
to this terrain. The sub-montane belt at the foothills of Siwalik along Himalaya is also 
very conspicuous, represented by high gradient piedmont zone, the Bhabhar and the 
adjoining relatively gently sloping Tarai belt. Bhabhar is a zone of unassorted sediments 
where fans and talus and disappearance of torrents traversing the southern scarp of the 
Siwalik are the common features, whereas, Tarai represents the zone of seepage. 
Bhabhar and Tarai also form source region of several streams like Gomati, Sukhata, 
Deoha etc. The rest of the area is a flat alluvial plain without any marked relief 
variations. The general slope of the rest of the plain is in south-easterly direction. 

2.2.2 Geology and Lithology of Ganga Plain 

The Ganga plain of U.P. comprises a thick pile of post Siwalik Quaternary 
alluvial sediments, overlying a northerly sloping basement (Sastri et al. 1971, 1973) 
comprising rocks of Bundelkhand Gneissic Complex, Delhi Group, Vindhyan Group and 
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Siwalik Group. The Peninsular Shield to the south exposes rocks of BundeUdhand 
Gneissic Complex, Delhi and Vindhyan Groups, whereas to the north of the plains are 
exposed rocks of Siwalik Group. The Quaternary alluvium of Ganga plain has been 
classified into Older alluvium coirprising Banda alluvium and Varanasi alluvium, and 
Newer alluvium consisting of Fan alluvium. Terrace alluvium and Recent alluvium 
(Bhartiya, 1996). 

Banda alluvium, comprising of numerous types of clays and red quartzo- 
feldspathic sands with translucent quartz, feldspar, chalcedony, chert, agate, jasper, 
fiugments of basic rocks etc. The next younger Varanasi alluvium, composed of oxidised 
silty clay and micaceous grey and brown sand, derived from Himalaya after “post Siwalik 
Epeirogeny” in Middle Pleistocene, the sedimentation continued up to Late Pleistocene. 
The Siwalik and a part of Banda alluvium which have been deposited prior to post 
Siwalik Epeirogeny form basement of Varanasi alluvium. 

Newer alluvium overlying BandaA^aranasi alluvium is represented by Fan 
alluvium near the foothill of Himalaya (Siwalik). The river valleys comprise with 
Terrace alluvium consisting of grey micaceous, laminated sand and thin beds of silty 
clay, and Recent alluvium comprising channel/bar deposits within the present bank limits 
of the respective rivers. Lacustrine deposits at places occurring over BandaA^aranasi 
alluvium also represent Newer alluvium. These Quaternary deposits have been classified 
broadly into three litho-stratigraphic units (Table 2-1) representing different cycle of 
sedimentation in basins of varying dimension and magnitude. 
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Table2.1 : Classification of Quaternary deposits of Gangetic basin (after Kumar 1996). 


Lithounit 

Distribution/Lithology 

NEWER 

ALLUVIUM 

Channel alluvium 

Active channel deposits of light grey, fine to 

medium grained micaceous sand and pebbles- 

cobbles (Himalayan provenance) and red Quartz- 
feldspathic medium to coarse grained sand 

(Peninsular provenance) along with drapes of silt 

and clays. 


Terrace alluvium 

Older Flood plain deposits (Terraces) composed of 

grey fine-grained micaceous sand with silt beds in 

Himalayan river. 


Alluvial Fan 

Deposits. 

Occur on the hill front and composed of sediments 

ranging from cobbles-pebbles to fine sand and silt. 

OLDER 

ALLUVIUM 

Disconfirmity 

Varanasi 

alluvium (VA) 

Polycyclic sequence of clay, silt and micaceous 

sand of Himalayan provenance with Kankar and 

ferruginous concretions 


Unconformity 

Banda alluvium 

(BA) 

(b) Chitrakoot Formation medium to course quartz 

feldspathic sand and silt clay of Peninsular 

provenance with minor kankar. 

(a) Variegated Clays 

PRE 

QUARTER- 

NARY 

Unconformity 

Vindhyan 

Supergroup/ 

Aravalli Group 

Bundelkhand Gneissic Complex in south and 

(Marwar Supergroup) 

Lower and Middle Siwalik groups in north 


I'he stratigraphic information obtained from deep wells in Ganga basin can be 


summarised as follows (Raju and Datta, 1960): 

1. The metamorphic basement is probably Aravalli schists and phyllites in the Badaun 
area. The nature of the basement may differ from one tectonic zone to another as 
could be visualised from the extension of peninsular Precambrian structural trends 
into Ganga basin, partly reflected in the aeromagnetic and seismic surveys. 








2. Overlying the basement, two major stages of sedimentation can be recognised. The 
lower stage is a typical platform association of the limestones-orthoquartzite-shale 
probably ranging from late Proterozoic to Palaeozoic. This association with well- 
defined folding probably extends as far north as the foothills zone showing a regional 
increase in the thickness. A regional unconformity representing a long period of 
hiatus in sedimentation (representing Late Palaeozoic, the entire Mesozoic and 
Paleogene) intervenes between these older sediments and continental Neogene sands 
and silts (Siwalik) representing the upper stage of sedimentation. This is particularly 
impressive all along the southern part of Ganga basin, whereas to the north, nearer the 
Himalayan foothill, the great thickness of sediments is more than 6 km. The neogene 
in Ganga basin represents gently dipping monoclinal starta towards the foothills 
without much structural deformation. 

3. The southern limit of Siwalik deposition progressively shifted southward with time 
along the southern margin of Ganga basin, the youngest Siwaliks or even alluvium 
rest directly over the Vindhyans or metamorphic basement. 

2.2.3 Geomorphology of the Area 

The area has been broadly classified into three major geomorphic units e.g.. 
Rocky surface, Upland plain and Flood plain (Figure 2.1). The further subdivision of 
these geomorphic units along with their elevations, landform elements and aerial 
distribution are given in Table 2.2. 



Figure 2.1 Geomorpho logical map of Ganga plain in U.P. (after Khan, 1996). 
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2.2.3.1 Rocky Surface 

Bundelkhand Gneissic complex surfece, Delhi surfece and Vindhyan smfece 
forming rocky surface of the south and constitutiag northern margin of the Pe ninsular 
Shield represent the oldest geomorphic unit of the area. Bunde lkhand Gneissic complex 
surfece, represent the subdued relief of the ancient most massif, comprises a few 
hummocks and inselbergs of granite generally 10 m to 80 m high above ground level and 
at times intruded by longitudinal wall like ridges of quartz reefe and dolerite dykes 
(Bisaria, 1996). Slopes of these ridges as well as the surro unding areas near foothill are 
regolith and colluvial materials and are mostly barren due to the absence or very little 
development of soil profile. 

2.2.3.2 Upland Plain 

It has been further classified into two sub-divisions e.g., Banda plain and Varanasi 
plain which represent oldest platform of Quaternary alluvium. 


Table 2.2: Geomorphic succession of Ganga plain, Uttar Pradesh (after Pandey, 1996). 


Geomorphic 

units 


Elevation 

above msl 

in meters 

Landform 

Elements 

Distribution 

Flood 

Plain 

Active Flood Plain 

40-300 

Point bars, 

channel bars, 

lateral bars 

Restricted within 

active channel 

domain of both 

Peninsular and Extra 

Penisular rivers. 


Other Flood Plain 

52-253 

Terraces, oxbow 

lakes, 

abandoned 

channels 

Confined to terraced 

valley of present day 

rivers. 


Alluvial Fan 

220-400 

Radial drainage, 

abandoned 

channels 

Near Himalayan 

foothill in patches. 










Varanasi Plain 

60-370 

Almost flat 

surface, tals, 

palaeo-channels 

oxbow lakes, 

sand ridges 

South-easterly 

sloping plain formed 

by palaeo-Ganga 

drainage system 

occupying most of 

the terrain north of 

Ganga Yamuna 

from Ghaziabad to 

Ballia 

Upland Plain 






Banda Plain 

90-160 

Undulating 

surface. 

Northerly sloping 

plain constituted by 

proto-drainage of 

peninsular rivers 

mainly confined to 

trans-Yamuna area 

between Agra and 

Mirzapur. 


Siwalik Structural 

Hills 


Sharp crested 

hills with steep 

obsequent and 

gentle dip slope 

From Saharanpur to 

Tanakpur all along 

the Himalayan 

foothills. 


Vindhyan Surface 


Plateau and flat 

topped isolated 

hillocks. 

Near Agra and 

Between Karvi and 

Mirzapur. 

RocIq' 

Surface 






Delhi Surface 


Ridges 

Between Delhi and 

Mathura 


BGC Surface 


Low lying dome 

shaped hillocks 

Between Jhansi and 

Karvi 


Banda Plain: It is an elongated tract covering parts of Bundelkhand and Baghelkhand 
tracks, the northern boundary of 'which is defined by Yamuna-Ganga axis and southern 
boundary by the Peninsular Shield. It mainly exhibits flat to moderately undulating 
topography with relief variations at micro level. Ravenous badlands stretching over 2 to 






3 km wide zones are developed on either sides of the main streams of this plain. The 
drainage pattern is coarse dendritic to sub-dendritic with a general north-easterly slope. 

Varanasi Plain: The vast alluvial terrain sprawling mainly north of Yamuna-Ganga axis 
up to foothill of Himalaya as well as narrow linear tracts lying marginally to the south of 
Yamuna-Ganga axis is termed as Varanasi plain. It represents flat to mildly undulating 
surface with a general south-easterly slope. Topographically the sub-mountain belt along 
foothill of Himalaya in the northern part of this plain representing Bhabhar and Tarai belt 
is a significant zone. Bhabhar belt, a piedmont, is a zone of unassorted sediments with 
fans and talus is a common feature. The Tarai belt represents a zone of seepage with 
shallow water table. The Bhabhar and Tarai belt also form the source region of several 
streams like Gomati, Sukhata and Deoha. The Bhabhar belt is a high gradient plain with 
slope varying from about 10 m/km in the northern part to about 0.40 m/km m the 
southern part aroimd Tarai Balt. Further, south-eastward the slope of the Varanasi plain 
goes on decrease which becomes 0. 1 m/km in the east. Ganga, Yamuna, Ghaghara and 
Ramganga are the major rivers draining the area which originate from the Himalaya 
while Chambal originating from Peninsular shield also forms a part of its drainage basin. 
The drainage pattern is coarse dendritic to sub-dendritic. Badland topography is observed 
along the either sides to main streams of the area. This surface comprises scores of fluvial 
signatures of paleo-drainage system in the form of tals, oxbow lakes, palaeo-channels, 
and paleo-levees, meander cutoff. 

Alluvial Fan: It forms part of the alluvial fan occurring at the foothills of Himalaya. It 
represents a piedmont zone formed due to coalescence of a number of fans. It mainly 
comprises coarser gravels, up to boulder size near its apex and sediments become finer 
(of sand-silt size) in the distal part. It shows steep slope and xmdulating topography near 
hill base but gentle slope and mildly undulating topography away from the hill base. 

2.2.3.3 Flood Plain 

The flood plain of present day rivers like Ganga, Yamuna, Ghaghara, Ramganga 
and Gomati comprises young geomorphic landforms, which include older flood plain and 
active flood plain of these rivers. 

Older Flood Plain: It represents terraced valley zones of the present day rivers and 
comprises both erosional and depositional terraces. Most of the rivers originating either 
from Himalaya or Peninsular shield, as well as alluvial upland have two terraces, one 



being erosional and the other depositional except Ghaghara where both the terraces are 
depositional. Erosional terrace also known as terrace of destruction (Fairbridge, 1968) is 
the older terrace among the two terraces and represent erosional bench of either Banda 
alluvium or Varanasi alluvium. This terrace occurring at relatively higher level is 
generally free from floods. It has moderately undulating topography and comprises 
abandoned channels, meander scrolls etc. Depositional terrace which represents, the 
lower level terrace, is occasionally affected by periodic floods. This terrace is made up 
of the sediments deposited by the streams to which it belongs and hence it is called 
terrace of construction. It exhibits mildly to moderately undulating topography with 
landform elements like abandoned channels, meander scrolls, cutoff meanders, crevasse 
splays, levees, and backswamps. In many rivers at places it has also developed a lower 
level sub-terrace component bearing resemblance in almost aU features and characters 
with the mother terrace. 

Active Flood Plain: This youngest geomorphic unit is restricted in the active zone of 
channel migration within the present bank lines of all the rivers. The least elevated 
surface represents active aggradational phase of the present day rivers comprising active 
landform elements like point bars, channel bars and lateral bars, a few of which are 
stabilised too. 

2.2.4 General Information about Study Area (Kanpur-Unnao) 

2.2.4.1 Introduction 

The Kanpur-Unnao region is a part of the central alluvial plain with its 
characteristic landforms. The area under study is located to the left side of Ganga river 
between Kanpur and Unnoa district. Kanpur, the eighth largest metropolis of India, is 
situated on the right bank of the Ganga river. It is the largest industrial town of the state 
of Uttar Pradesh. Unnao is situated on the left side of the Ganga river. The climate of the 
area is humid sub-tropical type. It is divided into three seasons: cold, hot, and monsoon. 
The highest temperature is about 42° C and the lowest is less than 10° C. The area 

receives annually 1 x 10®m^ of water/km* (Ansari, 1998). For the greater part of the year, 
high temperature prevails over the area and a large amount of water is lost by 
evaporation. 
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2.2.4.2 Subsurface Stratigraphy of Kanpur Area 

Bajpai (1989) has given the subsurface stratigraphy of the Marginal Gangatic 
plain of Uttar Pradesh using lithological information from borehole data of the Central 
Ground Water Board (CGWB) and other Groundwater Organisations. Figure 2.2 
represents a transverse (N-S) section from Kabrai in Hamirpur District to Unnao. This 
section shows basement at shallow depth (<100 m) up to Ghatampur in Kanpur District, 
it deepens suddenly near Panki (Kanpur District) and northwards, producing a highly 
variable alluvial thickness. 



Figure 2.2 Subsurface stratigraphy of alluvium in southern margin of 
Gangatic plain (redrawn after Bajpai et al, 1989). 

An important aspect of this subsurface stratigraphy is that in the northern part there is an 
important sandy (aquifer) zone between 270-370 m depth, which is made up of pink 




colour arkosic coarse sand and gravel (characteristic of all rivers coming from south). 
Below this zone are few horizons of fine-grained arkosic sands. These sandy horizons 
occurring in depths of less than 270 m are, however, of micaceous grey-sand type having 
Himalayan source. Between 270 to 300 m of depth sands show mixed sources. In the 
southern part, just above the shallow basement, a thick sequence of coarse-grained pink- 
coloured arkosic sand is followed on top by deposits of clay with kankar and fine-grained 
micaceous sand of Himalayan source. This subsurface stratigraphy demonstrates the 
following salient points: 

Initially, the rivers from Peninsular India were bringing fine arkosic sand and 
depositing on the slowly subsiding basement north of the present-day Ganga river. 
During this period the southern margin of the Gangetic plain was just south of Kanpur, 
and Ganga river was non-existent (Stage 1). 

A rejuvenation in relief, possibly due to thrust-fold loading in the Himalayan 
region caused an uplift of the Bundelkhand region (a peripheral upward) and rivers 
started bringing coarse-grained arkosic sand from the south (stage 2). The concurrent 
rejuvenation of relief in the Himalaya caused increased sediment supply from the 
Himalayan source and a southward shift in the depocentre. In Kanpur region, inter 
mixing of sediments from two distinct sources took place (stages 2 and 3) (Bajpai et al, 

1 989). Simultaneous activation of an E-W lineament caused gravity faulting and active 
subsidence of the area near Panki (Kanpur) in comparison to the area further south. 
Increased subsidence near Kanpur caused development of thick sand bodies (stage 4). 
Subsequently, the Himalayan source drainage and sediments, onlapped the southern 
source sediments, causing a southward shift in the southern drainage. The continued 
higher subsidence rate in the region north of the present-day Ganga river further 
enhanced the difference in the thickness of alluvium from south to north. For a long time 
Kanpur was an area where a major southern drainage met a major northern drainage, a 
situation analogous to Allahabad today. 

Figure 2.3 depicts the evolution of stratigraphy in Kanpur region sketched above. 
Though, no absolute dating is available, (Bajpai et al, 1989) proposed that the phase of 
coarse sand-gravel deposition near Kanpur corresponds to the Middle Pleistocene uplift 
in the Himalaya, when supposedly uplift of few thousand metres took place. 



1 () 



Figure 2.3 Schematic diagram showing evolution of southern margin of 
Gangatic alluvium (after Bajpai et al, 1 989). 





2.2.5 Tectonic Frame work of Ganga Basin 

A total of 15 boreholes touching the basement floor have been taken into 
consideration along section profile X-Y (Figure 2.5). The earlier studies have revealed 
that there are atleast ten NE-SW trending cross faults of significance which have not only 
affected the Neogene-Quatemary sediments but moved down to their basement floor 
(Figure 2.4) and thus have controlled the configuration of the basin. The studies have also 
shown that the entire stretch of the southern part of the Gangetic plain of U.P. between 
the AravaUi range in the west and the concealed Monghyr-Saharsa: ridge in east remained 
a positive land-mass till beginning of Late Pliocene (Figure 2.6). The faults 1, 4, 7 and 
10 are the oldest, possibly developed after Vindhyan sedimentation but prior to Banda 
alluvium sedimentation. Due to these faxxlts, the Vindhyan sediments are preserved in 
two grabens separated by a horst of Bundelkhvmd Gneissic Complex boxmded on either 
sides by faults 4 and 7. The strong tectonic movements during Late Pliocene, the fe,ult-10 
marking the eastern limit with Monghyr-Saharsa ridge, was reactivated, and a new feult-2 
developed in western part, along which the entire area lying between the two, down 
faulted to give rise to a basin for the deposition of Banda alluvium (Figure 2.6). The 
sedimentation in this basin continued till early Middle Pleistocene, when the entire area 
has been uplifted with reactivation of existing faults and development of new faults (5, 6, 
8 and 9). The entire Faizabad ridge acted as a graben but sub-divided into two grabens 
and a horst structure. Vindhyan graben east of Faizabad ridge bounded by feults 7 and 10, 
acted as a horst with a development of small graben limited on either side by faults 8 and 
9. The area west of Faizabad ridge has been relatively uplifted. 

The studies also pointed to pivotal nature of ftiults due to which the northern part 
of the Gangetic plain relatively subsided as compared to upliftment along southern 
margin of the basin. This has lead to the development of a new basin Indo-Gangetic basin 
for the deposition of the Varanasi alluvium commencing from Middle Pleistocene. This 
new basin laterally extended covering the entire Monghyr-Saharasa ridge in east and part 
of the Aravalli range in the west. 

The epiorogenic movements during Late Pleistocene times again changed the 
general scenario of this great alluvial basin. The Varanasi alluvium sedimentation 
culminated and was subjected to peneplanation, giving rise to Indo-Gangetic plain. 
During Holocene period, the Faizabad ridge fiirther subsided with upliftment of sub- 
block limited in fault 5 and 6. All the pre-existing faults except 9 have been reactivated 
and a new fault-3 developed in the western part of the basin. The area to the west of 



ridge uplifted in step like manner, along faults 1, 2, 3 and 4. The area to the east of 
Faizabad ridge lying between faults 7 and 8 has been uplifted, while that lying between 
faults 8 and 1 0 subsided giving rise to a large graben with upliftment of Monghyer- 
Sahrasa ridge. Fault-9 remained dormant during this period. 

The step like movement along faults 1 to 4 probably resulted in maximum 
upliftment of Aravalli range and development of south-easterly slope to the Gangetic 
plain controlling the flow trend of the rivers (Ghosh and Singh, 1988). The development 
of wide flood plain with Holocene Terrace alluvium in the Ganga river upstream of the 
fault-5 may be related to this step type of faulting where the western part of each block 
relatively subsided as compared to eastern one. The wide flood plain in the Ganga river 
down stream of fault-8 could be related to large graben developed west of monghyr- 
Sahrasa ridge. 



Figure 2.4 Sketch Quaternary geological map of Ganga basin showing locations 
of cross faults in the basement (after Khan et al, 1996). 
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2.3 GEOPHYSICAL INFORMATION ABOUT THE AREA 

The Indo-Gangatic basin is filled mostly with Late Tertiary and Quaternary 
sediments. A gradual thickening of the sedimentary layer from the Peninsular shield 
toward the Himalayan foothills was inferred by Oldham (1917) on the basis of gravity 
data. Recent, areomagnetic surveys, geophysical ground surveys and borehole data have 
indicated substantial variations in the thickness of the sedimentary overbtirden and 
imeven basement topography. Seismic surveys have been particularly useful in 
confirming the existence of buried ridges and basement upwards. 

2.3.1 Gravity Results 

A composite Bouguer anomaly map of Punjab and Ganga basins (Figure 2.7) has 
been prepared by using several sources (Choudhury, 1975 Balakrishna and Choudhury 
1979, Raiverman et al. 1983). Generally, the values are fairly reliable over the Gangetic 
plains. The contour interval is 10 m Gal. A few interesting observations from the gravity 
anomaly map are given below. 



Figure 2.7 Bouguer anomaly map over Ganga basin (after Verma, 1991) 



• Gravity contours have a WNW-ESE trend over the Punjab plains which becomes 
generally east- west over the Ganga basin. 

• Bouguer anomaly of the order of -50 m Gal has been found over the edge of the 
Peninsular shield and becomes increasingly negative towards the Himalayan foothills. 
Near the main boundary feult (MBF) in Punjab basin, the values are nearly -250 
mGal, while over the Gandak depression the minimum value is nearly -210 mGal. 
The gravity anomaly values become increasingly negative over the Lower Himalays, 
reaching -450 mGal over the Higher Himalayan region (Verma, 1991). 

• Gravity contours show large variations from Delhi towards Moradabad and Mathura 
towards Tihar, suggesting the presence of a ridge type structure underneath the 
alluvium or Upper Tertiary sediments^ 

• Gravity field does not show up as a prominent high over the Faizabad ridge. 

• Gravity field is much disturbed over Monghyr-Saharsa ridge, particularly in the area 
east of Patna. A part of this disturbance could be due to density variations within the 
crystalline basement. 

• A prominent gavity high is located over the Rajmahal hills. 

Profile Li L2L3L4 has been drawn nearly parallel to the axis of the Ganga basin, from the 
west to east. The basement depth (Karunakaran and Ranga Rao, 1979) and Bouguer 
anomaly have been shown along the profile as shown in Figure 2.8. This Figure shows 
that the basement depth is nearly 500-700 m below the surface in the area between Zira 
and Ujhani. Close to the Moradabad friult (Fi), the basement drops down from 750 to 
about 1000 m. Farther east, near Agra-Tilhar fault (F2), the basement shows a sharp 
change from 1000 m to 2500 m. The Agra-Tilhar fault appears to be a major friult, which 
has controlled the pattern of sedimentation in the western part of Ganga basin. The fault 
trends NE-SW from Agra towards Ujhani. Thereafter, it trends almost E-W towards 
Tilhar and beyond. The down-thrown block in this part is towards the south. East of this 
fault, the basement rises to nearly 1500 m along the fault (F3) and goes down to nearly 
2500 m over the western U.P. shelf. The gravity field shows a major low (-100 mGal 
value) over the western U.P. basin. East of the Lucknow fault (F4), the basement rises 
from 2500 m to about 800 m underneath the Faizabad ridge. The gravity field over the 
ridge shows no significant variations. East of Faizabad ridge the basement drops down to 
about 2000 m over the eastern U.P. shelf. The profile passes eastward over Patna fault 
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before entering the Monghyer-Saharsa ridge (MSR). The Bouguer anomaly, at a level of 
about -80 to -90 mGal over the eastern U.P. shelf, rises to about -20 mGal over the 
MSR. 

The major difference between the Faizabad ridge and the Ganga basin to its east 
and west lies in the thickness of the Upper Tertiary sediments. The Vindhyan as well as 
the Upper Tertiary sediments appear to be fairly thick underneath the western as well as 
eastern U.P. shelfs, compared to the area over the Faizabad ridge. The Vindhyan 
sediments, consisting mostly of quartzites, compact sandstones and phyllites (at places 
limestones), are of fairly high density (2.55 to 2.60 g/cm^) and do not contribute much to 
the gravity anomaly. The Siwalik sediments, being of low density (2.1 to 2.3 g/cm^), play 
a dominant role in controlling the gravity field. Although the gravity field over the 
Punjab and Ganga basins is influenced a great deal by the thickness of the Upper Tertiary 
sediments, this is by no means the only factor. The crustal thickening appears to 
contribute appreciably to the gravity field. 

PROFILE Li Lj Lj Li 



Figure 2.8 Profile L 1 L 2 L 3 L 4 nearly parallel to axis of Ganga basin, showing 
the major structures and nature of Bouguer anomaly. 

2.3.2 Seismic Studies 

Three regional seismic profiles have been carried out by Oil and Natural Gas 
Commission in the Ganga basin north of the exposed Vindhyan. These profiles are (i) 
Pilibhit-Hathras profile, (ii) Pilibhit-Bisalpur-Lucknow-Raebareli profile, and (hi) 
Nautanwa-Gorakhpur-Ghazipur profile. The seismic sections along these profiles have 
been extrapolated to the exposed Vindhyan outcrops towards south and south-west based 
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on geological investigations. The depth sections along the three profiles are shovra in 
Figure 2.9 (Rao 1973, and Karunakaran and Ranga Rao, 1979). In this section, the top of 
the Vindhyan formations below the Siwalik has been deiriarcated as an unconformity 
surface. It was not possible to demarcate the boundary between the Lower and Upper 
Vindhyan very clearly and, therefore, the two have been taken to represent together the 
total thickness of the Vindhyan sediments both in the outcrop axca covered by geological 
studies and in the Ganga Valley covered by seismic studies. Similarly, it has not been 
possible to demarcate clearly the boundaries between the Lower, Middle and Upper 
Siwalik and hence the three profiles have been taken to represent together the total 
thickness of the Siwalik sediments. On the basis of these three profiles, four maps have 
been prepared (a) Basement contour map of the Vindhyan basin, (b) Isopach contour map 
of the Vindhyan formations, (c) Contour map showing top of the Vindhyan, and (d) 
thickness of the Siwalik. Most of the geological sections, prepared by the ONGC, are 
based on the available seismic information. The seismic wave velocities obtained for 
different formations are given in Table 2.3. The velocity values have been converted into 
density using the well-known velocity versus density relationship obtained by Ludwig, 
and Nafe (Verma, 1991). 




Figure 2.9 Cross section based on seismic and geological data (after Verma, 1991) 



Table 2.3: Available seismic information for different formations 


from available sources (Verma, 1991). 


SI. 

No. 

Formation 

Location 

Depth 

(km) 

Velocity 

(m/sec.) 

Density 

(g/cm’) 

1. 

Upper Siwaliks 

Adampur-Hoshiarpur 

1.0 

2400-2600 

2.07-2.14 

2. 

Middle Siwaliks 

Adampur-Hoshiarpur 

2.0 to 2.2 

3000-3500 

2.2-2.3 

3. 

Lower Siwaliks 

Hoshiarpur 

2.0-3.0 

4500 

2.45 

4. 

Lower Siwaliks 

Janauri 

3.0 

4500 

2.45 

5. 

Lower Siwaliks 

Janauri-Gagret- 

- 

4500-4800 

2.45-2.50 



Dharmsala 




6. 

Lower Dhramsala 

Janauri 

5.0 

5400 

2.6 

7. 

Vindhyans 

Muzafl&rpur 

2.8 

5200 

2.5-2.6 

8 . 

Basement 

Janauri 

5.0 

6000 

2.7 

9. 

Basement 

Janauri-Balh area 

5.0 

6100 

2.72 


2.3.3 Aeromagnitic Survey in the Ganga Basin 

The aeromagnetic survey has been carried out over an area of more than 220,000 sq 
km of Ganga-Brahmaputra plains, and the depth to the basement has been estimated at 
different places distributed throughout the Ganga valley (Agocs, 1957). These results 
have been reviewed by Srivastava (1970). On the basis of qualitative and quantitative 
studies of aeromagnetic data the Ganga basin has been subdivided into several tectonic 
units separated by basement high (Figure 2.10). (Eremenko and Negi, 1968). These 
tectonics units are: 

• A shallow basin extending NW of Moradabad having a minimum depth of about 
5000 m of the basement. 

• A deep basin stretching from Moradabad to Lucknow (upto longitude of 82°). This 
maximum depth to the basement in this basin is around 6000 m. This basin is 
separated from unit (i) above by a basement feult trending NE-SW through 
Moradabad. The basement dips in the NE direction. 

• Basement uplift trending NE-SW in the Faizabad area may be the subsurfece 
extension of the Bundelkhand massif. 

• A deep basin stretching from NE of Faizabad up to Motihari. 

• A basin of intermediate depth east of Darbhanga and extending uto Pumea. 
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DETAILS OF IMAGE PROCESSING AND 
RESISTIVITY METHOD 


3.1 GENERAL 

This Chapter presents details of image processing and theory of data merging and 
brief discussion about edge extraction operators, which have been used in the present 
work. A brief discussion is provided on the resistivity method, which has been used to 
determine the subsurface resistivity variations. The present study involves following 
three major steps to achieve the objectives of the work: 

• Data merging technique 

• Edge detection 

• Resistivity Survey 

3.2 DATA MERGING 

The necessity of other sensors data for better resolution has attracted many 
scientists to develop data merging techniques. In this section data merging techniques are 
reviewed. 

3.2.1 Review of previous research 

Haydn et al. (1982) have used the first time the intensity hue saturation (IHS) 
procedure. They merged Landsat multi-spectral scanner (MSS) data with the return beam 
vidicon (RBV) data and Landsat MSS with heat capacity mapping mission (HCMM) 
data. The IHS method has also been used by Welch and Ehlers (1987) to merge thematic 
mapper (TM) and panchromatic (PAN) data of an urban area. They have concluded that 
the SPOT PAN data could be successfully used to enhance the interpretation of the 
Landsat TM mfirared bands. Thormodsgard and Fequay (1987) have used this method to 
merge spot multi-spectral and PAN data. Carper (1990) has suggested IHS method to 
merge SPOT, PAN and XS data. In this method, the intensity component was replaced by 
weighted average of PAN and XS3 ([2*Pan+XS3]/3) to improve the correlation between 
PAN and the intensity component. The study was conducted on three test sites. The first 


site contained both upland and lowland forest types and non-forested wetlands. The 
second test area consists of urban city and surrounding agricultural and forest areas. The 
third scene contained mainly bare fields with young com crops. Carper (1990) has 
evaluated his results by computing the correlation coefficients between original images 
and the images after data merging. He has concluded that IHS method is more suitable 
for areas covered by vegetation as the spectral values are better retained, whereas, the 
IHS method was quite useful for urban resource inventories. Jutz and Chorowicz (1993) 
have merged SPOT, PAN and Landsat TM data using IHS approach for detailed 
geological mapping. Harris and Murray (1990) have used IHS transform to merge radar 
imagery with other remotely sensed data. 

The high pass filter (HPF) method has been used by Chavez (1986) to merge TM 
data with digitised Panchromatic photograph collected as part of the national high altitude 
program data. The results were used to generate image products at scales up to 1:24,000. 
Chavez and Bo well (1988) have also used HPF to merge Landsat TM data with SPOT 
PAN data of agricultural, urban and geologic sites. They have compared the spectral 
information content of Landsat TM and SPOT data statistically usmg the correlation 
matrix and the principal component analysis, as weU as qualitatively through visual 
comparisons. They have concluded that the SPOT XS (multi-spectral) data for the most 
part duplicated the spectral information contained in the Landsat TM data, while the 
Landsat TM data have spectral information not contained in the SPOT XS data. 

Chavez et al. (1991) have performed a comparative study of IHS, PCA and HPF 
methods to merge Landsat, TM, SPOT and PAN data of an agricultural environment with 
a few small urban areas and a portion of the salt area. The results of the three methods 
have been compared in three ways: statistical, visual and graphic. Chavez et al. (1991) 
have concluded that the spectral characteristics of the data have been distorted in 
ascending order as follows: HPF, PCA and IHS. Pellemans et al. (1993) have used 
radiometric method to merge SPOT, PAN and XS images of an agricultural area and 
compared this method with IHS, statistical and spherical co-ordinate methods. They have 
concluded that the radiometric method preserved the absolute pbcel values (radiometric 
values) while the statistical method gives comparable results. The other two, IHS and 
spherical methods, provided more contrast in images which can be useful for the visual 
analysis. Munechika et al. (1993) have proposed a ratio method and have compared the 
method for merging Landsat, TM, SPOT and PAN data. They have computed the RMS 
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(root mean square) errors between the original and merged Landsat TM data and carried 
out a land cover classification, which recorded a 6% increase in the accuracy. 

Yocky (1996), Garguet-Duport et aL (1996), and Ranchin (1993) have used the 
multi-resolution analysis and wavelet transforms to introduce the spatial information into 
the spectral bands. Garguet-Duport et al. (1996) have con^jared this method with IHS 
and P+XS methods to merge SPOT, PAN and XS data for an application in landscape 
ecology. The study area consists of urban space, wetlands, and farmed zones with various 
yields and mining ponds. To interpret the performance of the method three criteria have 
been considered: geometrical improvement, preservation of the spectral characteristics 
and thematic analysis. It has been found that the wavelet method least distorted the 
spectral characteristics hence this method is well suited for the study of the vegetation. 
The visual comparison has indicated that the wavelet method generated results with good 
spatial resolution similar to IHS method or P+XS method. Yocky (1996) has used 
standard, additive and selective multi-resolution wavelet decomposition mergers to 
combine Landsat, TM, SPOT and PAN data and has compared these methods with IHS 
method. The comparison results have shown that the wavelet mergers provided greater 
flexibility and capability of preserving spectral values whereas, the IHS method found to 
change significantly. 

Vrabel (1996) has merged multi-spectral and panchromatic data from M-7 
airborne sensor. Study was conducted on areas of different combinations like residential 
and industrial, water treatment plant, storage tanks, agricultural, storage tanks, rail and 
desert roads, weapons test area and buildings. He has concentrated mainly on three 
points: validity of concept of effective GSD (ground sample distance), relative utility of 
band sharpening by different factors and comparison of IHS, HPF, P+XS and sparkle 
merging methods for various combinations dependent on the bands and the resolutions 
used. He has concluded that sparkle and P+XS methods give the best results in all the 
combinations tested. Wald et al. (1997) have used ARSIS, P+XS and duplication 
methods to merge SPOT, PAN and XS images of an urban area. They have proposed 
both, a formal approach and some criteria, to provide quantitative assessment of the 
merged images. Terrettaz (1998) has compared seven merging techniques: HIS, Carper 
(variant of IHS), HPF, P+XS, radiometric. Price and ARSIS (Ranchin, 1993) to merge 
SPOT, PAN and XS data in the context of an urban and sub-urban area. They have 
carried out detailed statistical and visual comparisons. It has been concluded that ARSIS 
and Price method (Susmitha, 1999) retain the spectral values but the visual improvement 



is not very good. The Carper (1990), radiometric and P+XS methods have also preserved 
the spectral values and achieved some degree of visual improvement but less than that of 
IHS and HPF methods. The IHS and HPF methods offer well-contrasted colour 
compositions that are very advantageous for the visual analysis. 

Following data merging techniques have been developed recently and details are 
given in recent literature. 

1 . High pass filter (HPF) 

2. Price algorithm I 

3. Intensity hue saturation (IHS) 

4. P+XS 

5. Principal component analysis (PCA) 

6. Radiometric integration 

7. Statistical method 

8. Spherical method 

9. Wavelet transform 

In the present study, four methods HPF, IHS, P+XS and Price have been used. 
The details of data merging techniques are given m the following sections. 

3.2.2 Data Merging Methods 

Data merging is combining of two data sets that are complimentary in nature. This 
is performed after the images get registered geometrically. This section explains in detail 
some of the available data merging techniques such as IHS, PCA, HPF, Price algorithm, 
radiometric, spherical and statistical techniques. The following description provides 
theoretical details of these methods. 

3.2.2.1 Intensity- Hue-Saturation (IHS) method 

The IHS method consists of transformation of low spatial resolution data set from 
the RGB to the IHS co-ordinate system and then the intensity component is replaced by 
the high-resolution data before carrying out the inverse transformation of the IHS to the 
RGB space. The use of this method rests upon the assumption that there is a close 
correlation between the high spatial resolution band and the intensity component. This 
assumption is valid if the high- and low-resolution bands cover in part at least the same 


spectral area. There are two models for representing colours: RGB and IHS models. The 
RGB co-ordinate system can be represented as a cube with red, green, and blue axes 
defining the X, Y and Z axes, respectively. The BBS model is hue-intensity-saturation 
model. The IHS co-ordinate system can be represented as a hexacone, where intensity 
represents total brightness and is represented by distance above the apex of the hexacone. 
Hue is the dominant wavelength of the colour and is represented by angular position 
around the top of the hexacone (0-360 degrees) and saturation is the purity of the coloiu: 
(i.e. the presence of white light in the image) and is represented by distance from the 
central, vertical axis of the hexacone. 

The advantages of IHS co-ordinate system over RGB system are summarised as 
follows: 

• The ability to vary each IHS component independently, without affectmg the others. 
For example, a contrast stretch can be applied to the intensity component of an image, 
and the hue and saturation of the pixels m the enhanced image will not be changed. 

• The information aspects of an image are presented using readily identifiable and 
quantifiable colour attributes that can be distinctly perceived. 

• In IHS system, numerical variations in the image data can be uniformly represented in 
easily perceived range of coloxirs. 

• The IHS approach may also be used to display spatially registered data of varying 
spatial resolution. For example, high-resolution data from one source may be 
displayed as the intensity component, and low-resolution data from another source 
may be displayed as the hue and saturation components. This advantage is the basis 
for using IHS transform as a merging technique. 

3.2.2. LI Steps Applied 

1. Transformation of three bands of the lower spatial resolution data set (slave images) 
from RGB to IHS co-ordinate system, the following important points should be kept 
in mind: 

• RGB values are converted to 0-1 before conversion to IHS coordinate. 

• IHS hue co-ordinate is rescaled from 0-360 degrees to 0-255 and saturation and 
intensity coordinates are converted from 0-1 to 0-255 before display. 

The Transformation from RGB to IHS system is carried out as described below: 



Since hue, saturation, and intensity are three properties used to describe colour, it seems 
logical that there be a corresponding colour model, IHS. When using the IHS colour 
space, one does not require knowing what percentage of blue or green to produce a 
colour. One simply adjusts the hue to get the colour one wishes. To change a deep red to 
pink, the saturation is adjusted. To make it darker or lighter, the intensity has to be altered 
to IHS. Many applications use the IHS colour model Machine vision uses IHS colour 
space in identifying the colour of different objects. Image processing applications-such as 
histogram operations, intensity transformations, and convolutions-operate on only an 
image intensity. These operations are performed much easier on an image in the IHS 
colour space. IHS is modelled with cylindrical coordinates. IHS, which uses the double 
cone model, looks like an inverted cone below another cone (Figure 3.1). 


white 



Figure 3.1 Double cone model of the IHS colour space 

The hue is represented as the angle 0 varying from 0° to 360°. The saturation 
corresponds to the radius, varying from 0 to 1. The intensity varies along the axis with 0 
being black and 1 being white. When S=0, the colour is a grey of intensity I. When S=l, 
the colour is on the boundary of top cone base. The greater the saturation, the farther the 
colotn is from white/grey/black (depending on the intensity). By adjusting hue the colour 
is varying from red at 0°, through green at 120°, blue at 240°, and back to red at 360°. 
When 1=0, the colour is black and therefore H is undefined, when S=0, the colour is grey 
scale and H is undefined in this case also. By adjusting I, a colour can be made darker or 



lighter. By maintaining S=1 and adjusting I, shades of that colour are created. The 
following formulas show how to convert from RGB space to IHS (Crane, 1997): 


I = j(R+G+B) 


5 = 1- 


(R + G + B) 


[min(i?,G,R)] 


(3.1) 

(3.2) 


H = cos ^ 


^(R-Gf+(R-B)(G-B) 


(3.3) 


If B is greater than G, then H=360°-H. 


2. In some cases contrast stretch is applied to high spatial resolution image so that it has 
approximately the same variance and mean as the intensity component image. 

3. Intensity component image is replaced by high spatial resolution image (master 
image). 


4. To convert from IHS to RGB, the process depends on which colour sector H lies in. 
For the RG sector (0" < H< 120°): 


b = 

r — 


-(1-5) 


1 + - 


5cos(/f) 
cos(60°-/f)J 


g==l-(r4-b) 

For the GB sector (120° ^^40°): 


H = H-\2Q° 


=1 Scos(H) 

3 cos(60" — H) 


(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 



'■=^(1--^) ' (3.9) 

b = \-ir + g) (3.10) 

For the BR sector (240°<H^60°): 


H = H-2A0° (3.11) 

b^lFu.. (3.12) 

g = |(l-S) ■ (3.13) 

r = l-(g + b) (3.14) 


The values r, g, and b are normalised values of R, G, and B, which are given below: 
R=3/r G=3/g 100B = 3/b 

Before display, R, G, B values are rescaled from 0-1 to 0-255. The algorithm has the 
limitation that only three bands can be merged at any one time. 


3.2.2.2 High Pass Filter (HPF) 

This HPF method merges the spatial information of the higher spatial resolution 
data set with spectral information of the higher spectral resolution data set. The high pass 
filter extracts the high frequency information from the image. The HPF method merges 
the spatial information of the higher spatial resolution data set with spectral information 
of the higher spectral resolution data set. In this method high pass filter is applied to the 
high spatial resolution data and the result is added pixel by pixel to the lower spatial 
resolution but higher spectral resolution data set. This allows the higher resolution image 
to be added equally to each of the lower resolution images using the simple pixel addition 
method without distorting the spectral balance. 

3.2.2.2.1 Steps Applied 

1 . High pass filter is applied on high spatial resolution satellite data. 

2. High pass filter can be applied directly on image to get high frequency 
information. Alternatively, one can first apply low pass filter such as mean filter 
and subtract this low frequency image from original image to get high frequency 
image. Procedure to extract high frequency information is explained below: 


a) If n x m be the size of the window of mean filter and iV x M be the size of the image. 



b) The window is placed m the top left comer of the image to be filtered and the average 
value of the elements in the input image that are covered by the window are 
calculated. This value is placed in the output image at the point corresponding to 
location of the central element of the window. 

c) In other words, the window can be thought of as a matrix with all elements equal to 1, 
the output from the convolution of the window and the image is the sum of the 
products of the corresponding window and the image elements divided by the number 
of elements in the window. 

d) If Wy, Gy, Oy represent elements in the mean filter, pixel values in the input image 
and pixel values in the output image respectively. Then 
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(3.15) 


where 

b= (n- 1)/2 
c= (m-l)/2 

n = number of rows in filter matrix 

m = number of columns m filter matrix 

r = central row in filter matrix (6+1) 

s = central column in filter matrix (c+1) 

i,j = image pixel underlying element (r, s) of the filter matrbc 


e) The size of the output image 0\s(N - 2b) x(M- 2c). 

f) Finally, the high frequency image is given by 

Hy=Gy-Oy. (3.16) 

3. Now add HPF results pixel by pixel to the lower spatial resolution but higher 
spectral resolution data set. 

Finally, we may get images with pixel values greater than 255 or less than 0. Pixel 
with grey value greater than 255 are assigned to 255 and pixels with grey value 
less than 0 are assigned to 0. 


4 . 



3,2.23 P+XS Method 

This method is based on the assumption that band P (Panchromatic) is the sum of 
the bands XSl and XS2. In P+XS method, the resultant merged XS images are given by 
equations (3.36) to (3.37): 


XS\+XS2 


(3.17) 


A5'1 + JSS2 


(3.18) 


PXS3 = XS3 


(3.19) 


3.2.2.4 Price a^orithm 

The objective of this merging method is to generate hybrid high-resolution multi- 
spectral images that attempt to preserve the radiometric characteristics of the original 
multi- spectral data. In the Price method (Price, 1987), each new band is the result of the 
combination of the original spectral band and an estimation of this band in terms of the 
high-resolution band. For the spectral bands closely correlated with the high resolution 
image, the estimation is obtained by linear regression, whereas, for weakly correlated 
bands a look-up-table is prepared to obtain the high resolution estimate. The principle 
behind the wavelet transform method is as follows: 

To pass from one level of resolution to a lower level, an image is decomposed 
into a context image, which represents the approximation of the image at high level 
resolution, and three detailed images representing the differences, in the horizontal, 
vertical and diagonal directions between the two levels. These three detailed images 
correspond to the spatial information of the high-resolution level. They are represented 
by the wavelet coefficients. If one knows this context image and the detailed images at a 
level of resolution, then one can reconstruct the high level image. For example, when 
IRS- 1C PAN and LISS-III images are merged, the wavelet transforms are used to extract 
the spatial structures between 6 and 24 meters in size from PAN and to introduce them 
into the LISS-III image. The LISS-III wavelet coefficients for resolution between 6 and 
24 meters are determined by models, which take into account the PAN wavelet 
coefficients and those of LISS-III at the other resolution levels. 

3.2. 2.4.1 Steps Applied 

Based on correlation between low-resolution and high-resolution data Price 
(1987) has proposed a two stage merger method. 



Stage 1: 

In this stage those low-resolution slave images are enhanced that are strongly 
correlated with the high-resolution master image and the procedure for this is as follows 
(Susmitha, 1999). 

1 . All pixels in the high-resolution master image that fall within a deiSned superpixel 
(area comprised by one pixel of slave image) are averaged and a new image 
(Mavg) is formed at the same resolution as the low-resolution slave image. 

2. The relation between pixel values of the slave images and master image can be 
described as linear. 

Si = aj . Move + b i (3.20) 

where 

S, Pixel value of the superpixel m the rth slave image. 

Move Average pixel value m the corresponding superpixel in the high resolution master 
iamge. 

a, , b, Least-squares regression coefficients. 

3. After solving for a,- and bi, a high resolution estimate of slave image can be 
formed as 

Sl=a,.Mhr + b, (3.21) 

where Mhr is the high resolution master image. 

4. High-resolution hybrid slave image is found out as 

Hybrid, = (3.22) 

‘^avg.t 

where 

Hybrid, Pixel value in the fth high resolution hybrid slave image 
.S',' High-resolution estimate of the ith slave image. 

S]^, Average of the pixel values m the S] image in the corresponding Si superpixel. 


5. 


Finally, merged images are these hybrid images. 



For an illustration, let m be the resolution of slave images and n = ml2 be the resolution 
of the master image where, 

(a) is slave image with resolution m 

(b) master image with resolution n 

(c) averaged master image with resolution w 

(d) high resolution estimate ( S] ) of the slave image with resolution n 

(e) averaged image of high resolution estimate (S"^, ) with resolution m 
(p) merged hybrid image (hybridj) with resolution n 

a, b, c, d, e, p, q, r, s, x, yi, y 2 , ys, y 4 , z represent pixel values. 

a + b + c + d , p + q + r + s 

e = and z = — — , 

4 4 

p = a^ia) + b, 

Similarly for q, r, and s\ where ai and b\ are regression coefficients, 
yx=-ip\y 2 =-iq\y^ =^-ir),y^ =-{s) 

Z Z Z ' z 

Stage 2: 

In this stage the weakly correlated slave images are enhanced. Here, relationship 
between the low-resolution slave images and the high-resolution master image can not be 
considered linear. In this stage a look-up-table is computed. This stage involves the 
following steps (Shaban, 1999). 

1 . All pixels in the high resolution master image that fall within a defined super pixel 
are averaged and a new image (A/avg) is formed at the same resolution as the low 
resolution slave image. 

2. To prepare a look-up-table, mean values of the pixel values in the weakly correlated 
slave images that fell in the same locations as the individual values in the averaged 
master image are computed. 

3. Referring to the high-resolution master image and look-up-table, place the mean 
values of the pixel values in the weakly correlated slave image as pixel values in high 



resolution estimate of the fth slave image i.e., Sj image in the same locations as the 
individual values in the high-resolution master image. 

4. Hybrid image of fth slave image is found out as 

Hybrid, = (3.23) 

'^avg.i 


where 

Hybridi^ixel value in the fth high resolution hybrid slave image 
5/ High-resolution estimate of the fth slave image. 

Average of the pixel values in the image in the corresponding S',- super pixel. 

5. Finally, the merged images are these hybrid images. 

(a) slave image with resolution m 

(b) master image with resolution n 

(c) averaged master image with resolution m 

(d) look-up-table to determine high resolution estimate 

(e) high resolution estimate {S ] ) of slave image with resolution n 

(f) averaged image of high resolution estimate (S^^, ) with resolution m 

(g) merged hybrid image (hybrid , ) with resolution n 

where, yi, y 2 , ys, y 4 are calculated in the same manner as that of stage 1 such that 

y, =~(.6),y, =-(4),y, =£(4).^. =-(6). 

5 z z z 

3.3 EXTRACTION OF EDGES 

One of the objectives of the present work is to emphasise the edges which are 
related to lithological contrast from an image using different operators and to understand 
the effect of data merging on extraction of edges especially on thresholding technique. 
So, it is necessary to have a brief idea on edge extraction procedure. The present section 
explains a brief idea about the edge operators, which has been used in this study. 


3.3.1 Introduction 

The edges of an im^e hold much of the information in that image. The edges tell 
where objects are, their shape and size and something about their texture. An edge is 
where the intensity of an image moves from a low value to a high value or vice versa. 
Information obtained from the images can be divided into two types, low frequency 
information and high frequency information. The slowly varying background pattern can 
be considered as low frequency information and the rapidly varying details can be treated 
as high frequency information. The filters that extract low frequency informations are 
called as low-pass filters and those extract high frequency informations are called high- 
pass filters. The high-pass filtered image can be used alone, particularly in the study of 
location and geographical distribution of edges. An edge is a discontinuity or sharp 
change in the greyscale value at a particular pixel point and it may have some 
interpretation in terms of geological structure or relief An edge in the digital image 
occurs at the boundary between two pixels when the respective grey level values of these 
adjacent pixels are significantly different. Edges are generally found because of changes 
in some physical and surface properties such as change in soil characteristics 
illumination, geometry (orientation, depth) and reflectance. An edge can be either 
physical or an intensity edge. A physical edge is one, which is actually present as a 
surface discontinuity whereas, an intensity edge is one, which is detected as sharp change 
in adjacent intensity values in an unage. Some intensity edges (like shadow caused by a 
strip of trees) do not correspond to an existing physical edge in nature. 

3.3.2 Edge detection 

The different approaches for spatial edge detectors can be categorised into linear, 
non-linear or best-fit methods (Peli and Malah, 1982). Linear methods involve edge 
detection based on differencing, such as first order derivatives evaluated in two directions 
(Rosenfeld and Kak, 1982) or in four directions (Pratt, 1977). Second order derivative 
operators like Laplacian and Laplacian of Gaussian (LOG) (Argyle, 1971; Macleod, 
1972) also come under linear approach. The non-linear approaches include Roberts and 
Sobel operators (Hall, 1979) or mask matching approaches such as Prewitt and Kirsch 
operators (Rosenfeld and Kak, 1982). The best fit methods involve fitting a simple 
function to the data and then using either the local gradient or deviations from it 
(Rosenfeld and Kak, 1982; Nalwa and Binford, 1986) or defining the edge as an 
intersection between different functions (Bums et al., 1986). These wide ranges of 



approaches have been summarised by Davis (1975), Rosenfeld and Kak (1982), Pratt 
(1977), and HaU (1979). Apart jfrom these, morphological operators are also available for 
the edge detection. Morphological operators are transformations of signals based on the 
theory of mathematical morphology (Cappellini et al., 1990; Haralick et al., 1987). Four 
morphological operations are erosion, dilation, opening and closing (Cappellini et al., 
1990). Levialdi (1983) has classified the edge operators based on their originating 
approach as: local, regional, global, heuristic, dynamic and relaxation. The following is 
the list of the various commonly used edge operators: 

• Linear operators 

• First order derivative 

• Second order derivative 
Laplacian 

Laplacian of Gaussian 

• Non-linear operators 
Roberts 

Sobel 

Prewitt 

Frei-Chen 

Kirsch 

• Morphological operators 

3.3.2. 1 Sobel Operator 

The Sobel edge detection masks look for edges in both the horizontal and vertical 
directions and then combine this information into a single metric. The masks are as 
follows: 


ROW MASK COLUMN MASK 


-_1 _2 -1 


0 

1 

1 

0 0 0 


-2 0 2 

_ 1 2 1 _ 


-1 0 1 


These masks are convolved with the image. At each pixel location we now have two 
numbers: Si, corresponding to the result from the row mask, and S 2 , from the column 



mask. We use these numbers to compute two matrices, the edge magnitude and the edge 
direction, which are defined as follows (Richards, 1993); 

EDGE MAGNITUDE 

4si+sl (3.24) 

EDGE DIRECTION 


tan"* 




*^2 


(3.25) 


The edge direction is perpendicular to the edge itself because the direction specified is the 
direction of the gradient, along which the grey levels are changing. 


3.3.2.2 Laplacian Operators 

The Laplacian operators described here are similar to the one used for pre- 
processing as described in section 2.2.3. The three Laplacian masks that follow represent 
different approximations of the Laplacian operator. Unlike compass masks, the Laplacian 
masks are rotationally symmetric, which means edges at all orientations contribute to the 
result. They are applied by selecting one mask and convolving it with the image. The 
sign of the result (positive or negative) fi-om two adjacent pbcel locations provides 
directional information, and tells us which side of the edge is brighter. 

LAPLACIAN MASKS 


1 

o 

1 

o 

1 

1 


■ 1 -2 1 ■ 


r-i -1 -11 

-1 

4 -1 


-2 4 -2 


1 

00 

! 

1 

o 

1 

o 

1 


1 -2 1 


-ij 


These masks differ firom the Laplacian type previously described in that the center 
coefficients have been decreased by one. With these masks we are trying to find edges 
and are not interested in the image itself An easy way to picture the difference is to 
consider the effect each mask would have when applied to an area of constant value. The 
preceding convolution masks would return a value of zero. If we increase the center 
coefficients by one, each would return the original grey level. Therefore, if we are only 
interested in edge information that is in the original image, the coefficients should sum to 
a number greater than zero. The larger this sum, the less the processedimage is changed 



from the original image. Consider an extreme example in which the center coefficient is 
very large compared with the other coefficients in the mask. The resulting pixel value 
will depend most heavily upon the current value, with only minimal contribution from the 
surrounding pixel values. 

3.3.2.3 High pass filter 

Sharpening an image is based on the high pass filter. A high pass filter will 
remove the low frequency components (like the image’s mean) and show only the details 
of the high frequency features. A common high pass filter ma.sk is 

999 999 999 

Another method of high pass filtering is to subtract a low pass filtered image from its 
original. This is known as unsharp masking. Another method for sharpening is high- 
boost filtering. A high-boost filtered image may be computed by increasing the intensity 
of the original and subtracting a low pass image. The equation is 

High-Boost = a Original - Lowpass 

when a=l, the result is a high pass filter. When a > 1, a fraction of the original is added 
back to the high pass result. This has the effect of restoring some of the low frequency 
components. The mask is: 


9 9 9 


where w = 9a-l 

9 9 9 


9 9 9 


Since high-pass filters enhance fine details in an image, they tend to amplify the noise in 
an image. This is seen in the output of high pass-based filters. 



3.3.3 Thresholding 

The final edge map is generated after thresholding the resultant image obtained by 
applying some edge operator. Thresholding refers to assignic^ all the pixel values below 
a certain level to ‘0’ and all the grey level values above the certain level to ‘255’ and that 
certain level is called as threshold value. Therefore, after thresholding a binary image is 
generated which contains only two pixel values ‘0’ and ‘255’. Pixels with grey value 255 
indicate an edge pixel. Thresholding will screen out all the unwanted variations in the 
edge image. Selection of a proper threshold value is a critical part. If we select a very 
high threshold value then it may cause the removal of most of the edges in the image. On 
the other side, if a low value is selected as threshold value, many spurious and unwanted 
edges will also be included in the final edge map. Therefore, an automated threshold 
selection is desirable than selecting a threshold value arbitrarily. 

Considerable literature is available for the selection of automatic Threshold values 
(Otsu, 1979). Majority of threshold detection methods use image histogram to optimise a 
certain criterion function. For a bimodal histogram, the threshold value may be taken at 
its valley bottom. In many cases, a single global threshold value may not give good 
results over the entire image. In order to overcome this problem, the image is divided into 
several parts and in each part threshold selection technique is applied. The final threshold 
value is nothing but the mode of the threshold values of all parts of the image. 

3.3.3.1 Otsu Thresholding Technique 

Image thresholding is typically used with edge detectors to emphasize the edges. 
This emphasizes the strong edges and de-emphasizes the weak edges. Thresholding can 
use one or two levels. Images can use one threshold value to set those above the 
threshold to the maximum pixel value and those below the threshold to 0. An upper and 
lower threshold can also be used to alter pbcel values. Values below the lowest threshold 
are set to 0. Those between the thresholds are left alone. Those above the upper threshold 
are set to 255. Thresholding technique is one that divides the image into two classes- 
background and objects. A survey on variety of thresholding techniques including both 
global and local thresholding was conducted by Sahoo et al. (1988). This section briefly 
explains a technique proposed by Otsu (1979) that has been used for the present work. 

Let V be the set total number of pixels, (x, y) be the spatial coordinates of the image, and 
G = (0, 1, , T-1 } be a set of positive integers representing grey levels. Then an image 



can be defined as a mapping ftinction f: NxN G. The grey level value of a pixel with 
coordinates (x, y) is denoted as/ (x, j;). Let teG be a threshold value and B = {bo, bj} be a 
pair of binary grey levels and bo, bjsG. The result of thresholding an image function at 
grey level t is a binary image function fr.NxN-^B, such that 


fi(x,y) = 



if f(x,y)<t 
if Ax,y)>t 


(3.26) 


In general, a thresholding method is one that determines the value t* of t based on a 
certain criteria. If t is determined solely jfrom the grey level of each pixel, then the 
thresholding method is point-dependent. If t* is determined from the local property (for 
example, local grey level distribution) in the neighbourhood of each pixel, then the 
thresholding technique is one that divides a given image into subimages and determines a 
threshold for each of these subimages. The Otsu thresholding method is described below 
(Otsu, 1979) 

Let the number of pixels with grey level i be n,. Then the total number of pixels in the 
given image is 


i-i 




(3.27) 


The probability of occurrence of grey level i is defined as 

P.=^ (3.28) 

i-i 

where p, > Oand = 1 

i=\ 

One has to divide the number of pixels of an image into two classes Co and C/ 
(background and objects) by a threshold at grey level t such that Q = {0, 1, 2, ...., t) and 
Cl = {t+l, t+2, ..., L-l}. Let cr^, cl, and be the within class variance, between- 

class variance, and total variance, respectively. An optimal threshold can be obtained by 
maximising one of the following conditions with respect to t 
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It can be noticed from Equations (3.70), (3,71) and (3.72) that cr^ and cr^ are functions 
of threshold level t, but are independent of t. It is also noticed that cr^ based on second 
order statistics (class variance) and cr\ is ba?ed on first order statistics (class means). 
Therefore, r\ is the simplest measure with respect to t and hence chosen as the criterion to 
evaluate optimal threshold at level t. Thus, optimal threshold t* is given below: 

t* = max Tf (3.39) 

3.4 RESISTIVITY METHOD 

3.4.1 General 

Of all the geophysical surveys, the electrical-resistivity method is employed for 
determining the electrical resistivity of the formation and also has found maximum 
applications in civil engineering problems. Electrical surveys are extensively used in 
civil engineering for determination the depth of bedrock, weathered zone thickness, 
groundwater occurrences etc. Electrical-resistivity surveys are easy to operate in the field 
with low expenditure. Electrical prospecting makes use of a variety of techniques, which 
include electrical resistivity, electromagnetic, induced polarization, telluric currents and 
magnetotelluric methods. The electrical-resistivity method is extensively used for 
groundwater, explorations. 

3.2.4 ELECTRICAL RESISTIVITY METHOD: 

Resistivity is defined as the resistance between opposite faces of a unit cube of the 

material. If a material of resistance ‘R’ has a cross-sectional area A and length L, the 
resistivity can be expressed as 



_RA 

P~ expressed in ohm-m 


(3.40) 


The electrical resistivity of the geological formation depends mostly on porosity, salinity 
and water content. All the rocks and soils transmit current due to the water content in 
void spaces. The resistivity of geological formation can be represented by the relation. 

Pf=Pw/FfFs (3.41) 

where Pf = resistivity of the formation 

Pw = resistivity of the water content in the formation ; 

Fw= Formation factor depending on porosity 

Rs = Resistivity depending on degree of saturation 
Hence the measurement of resistivity will be given as an indication of the presence of 
water and its salinity. In an electrical-resistivity survey, a known current I is sent into the 
ground through a pair of current electrodes Cj and C2 and the potential difference (V) 
created in the medium between a pair of potential electrodes Pi and P2 is measured. 
Hence, the resistivity of the formation based on Ohm’s law is stated by 

P-^Y (3.42) 

where K is termed geometric factor of the electrode arrangement 



<3.43) 


CiPi, P1C2, C1P2 and C2P2 are distances between the designated electrodes; P is the 

resistivity of the formation. However, the latter depends on the nature of the formation. 
If the formation is isotropic and homogeneous, the resistivity of the formation is 
measured. If the formation is anisotropic, apparent resistivity is measured. The apparent 
resistivity value depends on several variables, such as electrode spacings, geometry of the 
electrode array, true resistivity and geological factors. The effective depth of current 



penetration increases with an increase in electrode spacings. Various types of electrode 
configuration may be employed but Wenner and Schlumberger configurations are used 
more extensively for geotechnical and groundwater investigations. 

Cl P, P 2 C 2 


< a a a ► 


Figure 3.2 Wenner configuration 
(a) Wenner Configuration: 

The potential electrodes Pi and P 2 . are placed in a line with current electrodes Ci 
and C 2 . All are placed equi-distant from each other and disposed symmetrically with 
respect to a central point of investigation. A Wenner configuration is shown in Figure 
3.2. 


Pa = 2 n a— (3 44) 

where “a” is the distance between the adjacent electrodes, “I” is the applied current and 
“V” is the voltage difference between the potential electrodes. 


(b) Schlumberger Configuration: 

In this configuration all the four electrodes are placed in a line but the distance 
between the current electrodes Ci, C 2 is maintained equal to or more than five times the 
distance between the potential electrodes. A Schlumberger configuration is shown in 
Figure 3.4. 

The apparent resistivity is calculated by: 


p.=2n 


(L^-a^) AV 
4a I 


(3.45) 
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Figure 3.4 Schlumberger configuration 
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The depth of investigation in an isotropic and homogeneous formation is equal to half the 
distance between current electrodes. 

3.4.3 Field Procedure 

The resistivity surveys conducted in the field are mainly of two types: D.C. 
resistivity meter and A.C. resistivity meter. D.C. resistivity meter uses an external power 
supply, which can be increased to about 500 V and are capable of sending up to 2 
amperes. A.C. Instruments use a low-firequency alternating current source. 



Figure 3.5 Resistivity array 


Resistivity is measured by using four electrodes set in the ground. Current is 
applied to the ground through two of these electrodes and the resultant potential drop 
across the other two potential electrodes is noted as shown in Figure 3.5. Non-polarising 
copper sulphate porous pots are employed as potential electrodes to avoid various 
potentials developing between electrodes and electrolytes in the earth. A number of 
observations can be made by changing electrodes. The apparent resistivity is calculated 
from the known value of current (I) and the potential difference and the electrode 
spacings. The Wenner and Schlumberger methods are applied in the field. Field data is 
plotted on log-log graph sheet and the resultant field curve used for the interpretation. 


3.4.4 Resistivity Profiling 

There are two types of resistivity surveys, namely, profiling or lateral traversing 
and vertical electrical sounding (VES) or depth investigation. In a profiling survey, 
horizontal anomalies are delineated whereas in soundings vertical anomalies of the 
formation are investigated. 

(a) Resistivity profiling 

In this survey, electrode spacings are fixed; apparent resistivity values are 
determined at specified stations by shifting the whole electrode array along a particular 
profile or in a grid pattern. By measuring the apparent resistivity along the same profile 
with different electrode spacings and lateral as well as vertical changes in apparent 
resistivity can be obtained. 

Resistivity data collected during profiling is plotted against the corresponding 
stations along the profiles to depict an apparent-resistivity variations. If the readings are 
taken in grid pattern with one electrode spacings, the data may be utilised to prepare a 
contour map of apparent resistivity for the particular electrode spacings (Reddy, 1995 ) 

(b) Vertical electrical sounding 

In this method, keeping the place of observation constant, a set of apparent 
resistivity values are obtained successively for different electrode spacings. The value of 
apparent resistivity is plotted as a function of electrode spacings on log-log paper. The 
field curves may be grouped into two layers, three layers and four layers. The three-layer 

curves are represented as A, H, K and Q types. In A type P,<P2<P3, H type p,>P2<P3, K 

type p,<P2>p3, and Q type Pi>p2>P3. These curves are shown in (Figure 3 . 6 ). 

3.4.5 Interpretation 

The interpretation of resistivity data is carried out mainly in two stages. In the 
first stage the data is interpreted in terms of physical parameters, depth and resistivity. In 
the second stage interpreted results are correlated with geological knowledge to confirm 
the subsurface features. The popular method among interpretation techniques is curve 
matching. In this method the field curve is drawn on the same scale as the master curves 
to match with the master curves. In this study the interpretation has been carried out by 



using Resix standard resistivity package, which is based on the above idea. According to 
Reddy (1995), interpretation of three- layer curves is done as follows; 






Figure 3.6 showing the master curves for interpretation (after Reddy, 1995) 

1. Select the three- layer master curve sheet(s), which exhibits the same general shape as 
the field curve (A, K, H or Q). The master curves of each type are grouped into 

families having different resistivity ratio combinations (Pi:P 2 :P 3 )- The family of 
curves contains one curve for each h 2 / hi ratio. 

2. Select the three-layer master curve sheet(s), which exhibits the same general shape as 
the field curve (A, K, H or Q). The master curves of each type are grouped into 

families having different resistivity ratio combinations (Pi:P 2 :P 3 ). The family of 
curves contains one curve for each ha/ hi ratio. 
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3. Superimpose the transparent paper containing the field curve on the master curve 
sheet and by moving the sheets one with respect to the other (keeping areas parallel) 
to obtain a close match with one of the master curves. If an exact match is not 
obtained, interpolation between curves is made. 

4. Trace on the transparent paper the cross of the master curve together with the two 
resistivity marks Pi Pjand also the number which identifies the particular master curve 
(thickness ratio number). 

5. The thickness of the first layer will be equal to the cross abscissa on the field curve. 
The thickness of the second layer .will be equal to the thickness ratio number 
multiplied by the thickness of the first layer: depth to the top of the third layer will be 
equal to the sum of the thickness (first layer and second layer). 

6. The resistivity of the first layer will be equal to the cross coordinate on the field 
graph. The resistivity of the second layer will be equal to the coordinate to read on 
the field graph corresponding to the position of the resistivity mark. Similarly, the 
resistivity of the third layer corresponds to the p, resistivity mark. The interpretation 
procedure is shown in Figure 3.7. A perfect match is rarely obtained; interpolation 
between master curves gives satisfactory result. After comparing the field curves with 
standard curves formation parameters (thickness, resistivity, and number of layer) are 
obtained. Vertical electrical soundings are used for groundwater investigation and 
bedrock delineation. However, these surveys are also used for mineral exploration 
and toxic waste studies. 









Chapter IV 


RESULTS AND DISCUSSION 


4.1 GENERAL 

This chapter discusses the analysis of remote sensing data and resistivity data with 
special reference to deduce subsurface features. The remote sensing and resistivity data 
show very good correlation in studying the subsurface features. Efforts have been made 
to IRS- 1C LISS-in to improve to the visual interpretation and to understand the effect of 
data merging on thresholding technique for geological applications. 

4.2 METHODOLOGY 

For the mapping of surface features, we have used remote sensing data and for 
subsurface mapping resistivity soundings have been carried out. The approach followed 
in the present study is shown in Figure 4.1 : 

4.2.1 Image processing: 

The following steps have been carried out through this study: 

1 . Generating false colour composite (FCC) using blue, green, and red bands. 

2. Generating Normalized Vegetation Index (NDVI), and reclassification of image 
based on NDVI values. 

3 . Extraction of geomorphological features from stretching FCC. 

4. Applying edge detection operators (Sobel, Laplacian, High pass, and Otsu 
thresholding technique) on green and red bands. 

5. Registering of multi-spectral images to Pan characteristics using Idrisi package. 

6. Merging of the registered images using the HPF and Price algorithm, IHS, and 
P+XS methods. 

7. Final image showing fault location. 

4.2.2 Resistivity soundings: Resistivity soundings at 14 locations have been carried out. 
Following steps have been followed to deduce subsurface resistivity information: 

1. We have calculated the apparent resistivity using P ; while: 



Where C and P are positions of current and potential electrodes, respectively 
discussed in Chapter IE. 

2. Analysis of resistivity data using Resix package. 



Figure 4.1 shows the methodology followed in present study 




















4.3 Digital image processing: 

Detailed processing of digital remote sensing data have been carried out to deduce 
surface features. 

4.3.1 False Colour Composite (FCC): 

Using multi-spectral bands from IRS- 1C LISS-III, MiR , green, and red bands, the 
FCC image has been generated (Figure 4.2). IRS-IC LISS-III false colour composite 
image shows a channel migration pattern through neck cut-off spread over an area about 
150 Km^ along the Ganga river in the (Bithoor) Kanpur region. FCC shows different 
areas representing agriculture, urban areas, and geomorpho logical features like prominent 
meander cut-off. Ox-bow lake, and successive of paleo-channels. 



Figure 4.2 FCC false colour composite of the study area 
4,3.2 Geomorphological studies: 

From stretching of FCC (Figure 4.3), we have generated geomorphological map 
(Figure 4.4) of the study area. Geomorphological studies have been carried out, which 
indicate that the Ganga river has reached its present location through 6 steps, in the first 
step the Ganga river was flowing at the eastern side of the Unnao city. Then started its 
migration in western direction to flow at west of Unnao, this has happened in second 



step. In the third stage the northern main channel has shifted to the left and has created a 
loop channel. In the fourth step, the Ganga river has migrated towards eastern direction 
may be due to subsurface movement and started its westward movement which lead to 
meander cut-off in the fifth step. At the last step, river changes its main channel to 
present course. Inside meander cut-off, the image clearly show the presence of two 
distinct geomorphic zones, in the western zone older formation of flood plain 
characterised by modified or overlapping paleo-channels. The eastern part characterised 
by development of meander scars, sloughs, and lakes. In this zone younger deposition is 
found compared to the western part. It is clearly seen from the image that the 
geomorpho logical features like paleo-channels show discontinuity and change their 
direction. The change in direction of these paleo-channels reflects the upward and 
horizontal movements of the subsurface. 






c=2a Paleo channel 
^ Water body 
^ ^ Loop boundary 



Figure 4.4 Geomorpho logical map produced from FCC. 


4.3.3 NDVI Normalized Vegetation Index: 

Using NIR and Red bands, we have computed NDVI using the following formula: 

_ Band4(NIR) - Band3(Red) 

Bands + Band4 

NDVI image is shown in Figure 4.5. This NDVI image is further classified in seven 
classes representing water body, very poor, poor, moderately, high, very high vegetated 
areas, and tree with 2 meter height. Classified image is shown in (Figure 4.6). 


Table 4.1 NDVI values for different classes 



NDVI Values range 

Vegetation condition 

Class I 

-0.362 - 0 

Water 

Class II 

0-0.196 

Very poor 

Class III 

0.196-0.296 

Poor 

Class IV 

0.296-0.396 

Moderately 

Class V 

0.396-0.496 

High 

Class VI 

0.496-0.596 

Very high 

Class VII 

0.596-0.696 

Tree with 2 meter height 


The normalised difference vegetation index (NDVI) values show that the study 
area can be classified in two distinct vegetated zones: western zone is highly vegetated 
than the eastern part, which shows negative values for NDVI due the presence of water, 
and poorly vegetated. The highly vegetated areas are related to paleo-channels. 










4.3.4 Edge detection: 

One of the main objectives of the present work is to know the effect of merging 
on edge detection. To compare the results, edge extraction procedure involves two major 
steps on both the original and merged data. Band 2, which is found to suited more to 
extract geological features compared to Band 3, which shows more number of linear 
features. Band 2 is used for extracting edges of lithological contrast in both cases before 
merging and after merging. These images are shown in Plates 4.1 to 4.6 for IRS-1 C 
Bands 2 and 3. For the extraction of the edges following steps have been followed: 

(a) Edge detection 

1 . Use of Laplacian, second order derivative operator on Band 2 image to compare with 
Sobel operator image. The output images are shown in Plate 4.1 which shows that 
Laplacian operator is unable to give good details in studying edges which are related 
to lithological contrast as given by Sobel operator. 

2. The Sobel non-linear operator, which is considered as one of the best available 
operators, has been used to extract the edges from the image. The output images have 
been generated (Plate 4.2) using absolute gradient values. 

3. Use of High passes filter to remove the low frequency components. The high 
frequency features are shown in Plate 4.3. 

(b) Thresholding 

The Otsu technique has been used to find the optimal threshold value. The image 
is divided into 1 6 different equal parts. Otsu technique is applied to Bands 2 and 3 before 
merging and after merging. The global threshold value obtained as the mode of threshold 
values for all parts in the image. Finally, the global threshold values one for each band m 
two cases have been obtained. Using this threshold value the image has been converted 
into binary form by assigning the pixels value less than the threshold value to ‘0’ and 
pixels with grey value greater than the threshold value to ‘255’. Otsu thresholding 
technique has been applied on Bands 1 and 2 before and after applying merging 
techniques using only HPF, and Price methods. Otsu thresholding technique has given the 
ability to map the fault locations and to recognize the high and low land areas. Different 
values of thresholding might be useful in numerous geological applications. Table 4.2 
shows the optimal thresholding values for the two cases. 


Table 4.2 asu threshold value for IRS-IC Band 1 and 
Band 2 using 16 image parts 


Multi-spectral 

bands 

Before merging 

After merging 

Thresholding values 

HPF 

Price 

Thresholding values 

Thresholding values 

Band 2 

37 

16 

18 

Bands 

48 

18 

26 


From the optimal threshold values, it is clearly seen that merging technique has 
filtered the noise due to which more details have appeared in the output images. It has 
been found that the use of Otsu thresholding method to a merged data has not given better 
results compared to the original data. Applying Otsu technique on Band 2 we have foxxnd 
better geological features and clear image with more information about the area 
compared to those of Band 3 (Plates 4.3 to 4.6). 

4.3.5 Merging Technique 

After registration, IRS-IC LISS-III multi-spectral images to the PAN data, 
different merging algorithms i.e. IHS, P+XS, HPF and Price have been applied. The main 
aim of the merging techniques is to introduce the spatial information of the high spatial 
resolution PAN image to the low spatial resolution images and IRS-IC LISS-III data 
without distorting the spectral characteristics of the higher spectral resolution bands and 
IRS-IC LISS-II data. In the present study, the merging of IRS-IC LISS-in with PAN 
data has been carried out. 

The following discussions briefly explain different input and output data used in 
various merging algorithms. 

(a) IHS 

In this algorithm, only three images can be merged at a time. The registered multi- 
spectral Bands 1 to 3 are merged with PAN data. The registered low-resolution images 
are transformed from RGB space to IHS space. Before taking inverse transform, the 
intensity component is replaced with the PAN data. The main assumption for replacing 
the intensity component with the high-resolution data is that the two images are 
approximately equal to each other spectrally. The merged images are shown in Plates 4.7 


and 4.8. 
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Laplacian edge operator applied on Band 3 


Plate 4.1 






Sobel edge operator applied on Band 2 



Plate 4.2 





Otsu thresholding applied on Band 2 before merging 
using optimal value equal to 40 


Plate 4.3 



Otsu thresholding applied on Band 2 before merging 
using optimal value equal to 37 



Otsu thresholding applied on Band 3 before merging 


using optimal value equal to 48 


Plate 4.4 




,\.r. '■ c i,- i ’ p- ^ 

^4. <•*' .•' V ' ' ' r"-''*'' ■* •"■ ' ^■''■^/'U.9 


Otsu thresholding applied on merged Band 2 using HPF 
merging method using optimal value equal to 1 6 
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Otsu thresholding applied on Band 3 using HPF technique 
and optimal value equal to 1 8 
Plate 4.5 



Otsu thresholding applied on merged Band 2 using 
Price method and optimal value equal to 18 



Otsu thresholding applied on merged Band 3 using 
Price method and optimal value equal to 26 
Plate 4.6 
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(b) P+XS 

In this method two images of Band 1 and Band 2 have been merged. Band 3 is 
same as before merging. This method is based on the assumption that Band P 
(Panchromatic) is the sum of the Bands 1 and 2. First multi-spectral images of Band 1 
and Band 2 have been registered to Pan band. The output images are shown in Plates 4.8 
and 4.9. Using equations 3.17 to 3.18 the resultant merged images have been generated. 

(b) HPF 

HPF (high pass filter) is applied to high spatial resolution data, and mean filter is 
applied on low resolution data and the results are added pixel by pixel. The size of the 
mean filter applied on the PAN data is approximately equal to the ratio of resolution of 
different satellite data used for merging. Hence, mean filters of size 3 x 3 is applied, 
respectively on PAN data and the resultant image is subtracted from the original image to 
extract the high frequency image. This is added pkel by pixel to each low-resolution 
slave image. The merged images are shown in Plates 4.10 and 4.1 1. 

(d) Price 

This method uses the original slave images and registered or without registered 
images. The Price algorithm is based on the correlation between high-resolution and low- 
resolution images. The correlation is determined between these images, PAN data is 
either averaged or registered to multi-spectral bands. The correlation values are given in 
Table 4.3. As there is a difference in the date of acquisition of the PAN, and multi- 
spectral images of IRS- 1C LISS-III, experiments have been carried out using only PAN 
image registered to multi-spectral images using only nearest neighbour resampling 
techniques. There are two steps for merging with the Price algorithm based on the 
correlation coefficients. In step-1, the highly correlated slave images are used for 
merging. The step-2 is applied for weakly correlated slave images. 

The first step of Price algorithm is used to merge Bands 1 and 2. The second is 
used to merge Band 3. In Table 4.3 the correlation coefficients values for Bands 1 and 2 
are small because of different dates for Pan and multi-spectral images, while in Band 3 it 
is normal. Table 4.3 also shows that the use of registered Pan image gives better 
correlation coefficients values compared to the use of averaged Pan image. The merged 
images are shown in Plates 4.1 1 and 4.12. 
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Table 4.3 Correlation coefificients between PAN and Low Resolution images 


Low resolution 
images 

PAN^ 

PAN,,g(NN) 

Band 1 

0.1196 

0.1886 

Band! 

0.1286 

0.1973 

Bands 

0.0434 

0.1263 


PANavg - averaged PAN data 
registered PAN data 

4.3.5. 1 Results from Data Merging 

Using four methods (HPF, Price, IHS, and P+XS) the visual comparison have 
been made between original and merged data. Using the above mentioned methods give 
the following important points: 

1 All methods are found to improve the resolution as well as features present in the 
image. 

2 Visually HPF and Price techniques have given better results compared to other 
methods. 

3 Only HPF and price method have removed cloud cover from output merged 
images. 

The root mean square errors (RMS) are the measure of the accuracy of the Ground 
Control Points (GCP) selected for registering the images. The RMS error is expressed in 
the input image units in the form of number of pixels which should be less than 3 pixels. 
The results obtained from the experimental studies are presented in Table 4.4. No 
significant effect of RMS values to merging results are seen from Table 4.4 because it is 
obtained equally for different merging methods and with considering different bands, 
except for Price method. 


Table 4.4 Average RMS errors of GCP 



RMS errors used for 

RMS used 

Bands 

HPF,1HS and P+XS 

for Price 

Band 1 

0.592 

0.45 

Band 2 

0.592 

0.45 

Band 3 

0.592 

0.45 
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Statistical and numerical comparisons have been made to evaluate the effect of 
different data merging techmques on distortion of the spectral characteristics of higher 
spectral resolution data. The mean and standard deviations have been calculated before 
and after merging for different merging technique as well as for different bands. The 
mean and standard deviations values are given in Table 4.5. 


Table 4.5 Mean and Standard deviation values of original and merged IRS-1 C data. 


Band 

Estimates 

Original 

Data 

Merging methods 

Nearest neighbour 

HPF 

Price 

ms 

P+XS 

Band 1 

Mean 

104.386 

100.94 

103.7 

136.34 

158.6 

S.D. 

13.173 

9.787 

11.22 

27.9 

BMiHI 

Band! 

Mean 

67.856 

65.82 

87 

132.4 

IREIiH 

S.D. 


10.894 

9.482 

29.11 

20.87 

Band 3 

Mean \ 

92.660 




- 

S.D. 

11.593 

10.667 

9.898 

23.98 



An absolute difference image between the merged and the original image is 
computed. The mean and standard deviation of the difference image indicates the 
distortions due to merging. Higher the value of these factors, more the distortions due to 
merging technique used. The number of zero pixels present in these (difference) images 
are also determined. Higher the number of zero pixels, lesser the distortion in the spectral 
values. The number of zero pixels has been calculated from difference images between 
original and merged data. Table 4.6 shows how the merged data differ from the original 
data numerically. 


Table 4.6 Number of zero pixel values of merged IRS- 1C data. 


Band 

Merging methods 

Nearest neighbor 

HPF 

Price 

HIS 

P+XS 

No. of ‘0 'pixels value % 
for Band 1 

40% 

10% 

2% 

0% 

No. of *0 'pixels value % 
for Band 2 

41% 

6% 

0% 

Less 
than 1% 

— 

1 

No. of *0 'pixels value % 
for Band 3 

38% 

8% 

1 

Less 
than 1% 

- 


From Mean, Standard deviation, and Number of ‘0’ pixel values, it has been 
found that the HPF method gives the best results with less distortion to spectral 
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characteristics. The number of ‘0’ pixels obtained for different bands is around 40%, 
while from Price techmque it is found around 9% only, this may be related because of 
the use of different date for Pan, and multi-spectral Bands. Generally, the other methods 
IHS, and P+XS methods give number of ‘0’ pixels less than 1%. The mean and standard 
deviation values of the merged data from the HPF and Price algorithm are numerically 
very close to that of the original data. The mean and standard deviation values of data 
resulting from the IHS and P+XS methods are much deviated from the original values. 
Based on the numerical closeness of the mean and standard deviation values of the 
merged data to that of the original data, different merging techniques can be arranged in 
the descending order of preference as HPF, Price, IHS and P+XS. 

Following points have been adopted while applying merging methods: 

1 . Both HPF and Price algorithm have given very good results but it is better to use 
HPF method if there is difference m time for data acquisition between 
Panchromatic and multi-spectral bands. 

2. For regional geological investigation it is not necessary to use merging technique 
but for study the features from a specific region one must use merging techniques. 

3. Use of Price and HPF methods remove the cloud cover from the merged data. 

4. Use of merging techniques reduce the noise and give more details of an image 
which are shown from optimal Otsu thresholding values. 

5. The applied merging techniques have not given any improvement for geological 
edge extraction the reason is related to scale and its effect on pixels from merged 
data. 

4.3.6 Results from remote sensing data: 

From the detailed remote sensing data analysis following results have been 
obtained. 

1. I'he geomorphological features show the inside of meander cut-off image clearly 
showing the presence of two distinct geomorphic zones: the western zone is 
covered by older flood plain characterised by modified or overlapping paleo- 
charmels. The presence of two distinguished areas may be related to major 


tectonic events. 



2. The vegetation index (NDVI) shows two different areas, western zone with 
highly vegetation compared to other parts, which show negative values for NDVI 
due the presence of water, and poor vegetation. 

3. Sobel edge operator also shows two distinct zones: western part with more edges 
compared to that of eastern part. 

4. Otsu thresholding technique is found to show clearly the presence of lineaments 
representing faults in the area and also gives information about the presence of 
low land area in the eastern part and uplifted area in the western part. 

5. Band 2 image clearly shows edge enhancement compared to Band 3 image, and 
also even after applying merging techniques. The edge operators on original data 
found to give better results compared to merged data. 

6. Merging technique helps to understand the remaining edge of the channel 
migration due to presence of road not related to the presence of fault. 

7. For studying subsurface structure, it might be usefiil to study topography 
variations through generating Digital Elevation Model (DEM) and draped image. 
In the present study due to unavailability of complete topography data it was not 
possible to use this approach and its relationship with subsurface structure. 

8. The subsurface and surface features from remote sensing data are shown in Figure 
4.7. 











Band 1 after merging using High Pass Filter (HPF) method 



Band 2 after merging using (HPF) method 


Plate 4. 10 





Band 1 after merging using Price method 


Plate 4. 1 1 




Band 3 after merging using Price method 


Plate 4.12 
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4.4 Geophysical Data Analysis 
4.4.1 Resistivity method 

The electrical resistivity of rock formation limits the amount of current passing 
through the formation. It is the resistance in ohms between opposite faces of a unit cube' 
of the material. The resistivity of rock formation varies over a wide range, depending on 
the material, its density, porosity, pore size and shape, water content and quality, and 
temperature. There are no fixed limits for resistivity of various rocks: igneous and 
metamorphic rocks vary in the range lO’to 10* ohm-m; sedimentary and unconsolidated 
rocks, 10" to 10' ohm-m. Figure 4.8 provides a representative guide to the electrical 
resistivity ranges of various sedimentary rocks (Todd, 1980). In relatively porous 
formations, the resistivity is controlled more by water content and quality within the 
formation compared to those by the rock matrix. The resistivity decreases with the degree 
of saturation and salinity of the ground water. The clay minerals conduct electric current 
through their matrix; therefore, clayey formations tend to display lower resistivity than 
compared to permeable alluvial aquifers. 
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Figure 4.8 A representative guide to resistivity of various 
sedimentary rocks (after Todd, 1980). 




4.4.1.2 Field work 


Resistivity survey has been carried out in Kanpur-Unnao region near Bithoor at 
the left bank of Ganga river. Using OYO resistivity meter 14-resistivity soundings along 
the meander cutoff have been taken. These sounding points cover the entire area of 
channel migration pattern. The sounding locations are shown in Figure 4.9. 



Figure 4.9 Location of resistivity soundings on FCC image. 


The readings have been taken using Schlumberger configuration, which is 
described in chapter El. The field procedures (Figure 3.5) have been discussed in chapter 
IE. At some locations (Dostinagar, DuUakhera, and Garha) it was difficult to have current 
electrode spacings more than 100 to 150 m due to the presence of vegetation cover. For 
most of locations current electrode spacings extended up to 600 m. In Figure 4.10 the 
OYO resistivity meter and its set up are shown. 


Figxire 4.10 Resistivity set up in the field. 


4.4.1.3 Results and discussion 

After field work, the detailed processing of data have been carried out. Th 
apparent resistivity have been plotted with half current electrode spacings (AB/2) usin 
semi-logarithmic sheet. Resix package has been used to interpret these curves. Fror 
Resix package, all the parameters including number of layers, true resistivity, and th 
thickness of each layers have been obtained. The variations of apparent resistivity alonj 
profiles I, H, ni, IV and V are shown in Figures 4.11, 4.12, 4.13, 4.14 and 4.15. Profil 
(I) represents the aparent resistivity variations from south to north direction of study are 
along Dostinagar, Garha, Kima, and Firou 2 pur (Figure 4.11). Profile (II) represents fin 
aparent resistivity variations from south to north direction along Dullakhera, Parmani 
Diwankhera, and Isunyian (Figure 4.12). Profile (III) represents the aparent resistivib 
variations from south to north direction along Utaria, Mahipatikhera, and Powa (Figur( 
4.13). Profile (IV) represents the aparent resistivity variations from south to nortl 
direction along Pyarepure, Bhadeona, and Mubarkpur (Figure 4.14). Profile (V 
represents 3 -dimensional variations of aparent resistivity from east to west directioi 
along Firuzpur, Isunyian, Mubarkpur and Mahipatikhera (Figure 4.15). The contour maj 
has been generated from apparent resistivity values using 30 ohm-m contour intervals 
this map is shown in Figure 4.16. By studying the profiles I, n, IE, IV, V, and contou 
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map, low apparent resistivity is clearly seen at locations Dostinagar, Dlluakhera, Garha, 
Pyarepure, Utaria, Parmani, Kima, Powa, and Diwankhera. While at all the other points 
high apparent resistivity values are found. The highest resistivity values up to 500 ohm-m 
are found at locations Firu 2 pur, Isunyian, Mubarkpur and Mahipatikhera (Figure 4.16). 
Only at Dostinagar apparent resistivity shows gradually increase which may be related to 
the location of Dostinagar which is lying far away to the channel migration pattern 
(Figure 4.9). An interpretation of resistivity data using Resix package shows a very good 
contrast between the eastern and western parts of the meander cut-off. The layered model 
estimated from Resbc package for each point and give very good fit with the field curves 
(Figures 4.17 to 4.30). The data estimated shows a large contrast in resistivity values 
between two sides, while they are varying between 2-60 ohm-m at locations Dostinagsr, 
Dullakhera, Garha, Pyarepure, Utaria, Kirna, Powa, and Diwankhera. The low resistivity 
values indicate the presence of clayey, shale and fined sand formations. The amount of 
sand and the presence of hard shale may give rise to the resistivity up to 60 ohm-m. 




Figure 4.1 1 - Profile I (S-N): Apparent resistivity variations along 
Dostinagar, Garha, Kirna, Firouzpur. 



Figure 4.12 -Profile II Apparent resistivity variations along 
Dullakhera, Parmani, Diwankhera, Isunyian. 



Figure 4.13 - Profile III (S-N): Apparent resistivity 
variations along Utaria, Mahipatikhera, Powa 










Figure 4.14 Profile IV (S-N): Apparent resistivity variations along 
Pyarepure, Bhadeona, Mubarkpur. 



Figure 4.15 - 3-diniension Apparent Resistivity variations 
along Profile V Firu 2 pur, Isimyian, Mubarkpur and Mahaptikhera. 
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Figure 4.16 shows Apparent Resistivity contours in the study area. 
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The data collected from the eastern part was easy to interpret since the fitting errors 
between master and field curves was very small compared to locations Firu 2 pur, 
Isunyian, Mubarkpur, and Mhaiptkhera in the western part. In the eastern part the 
resistivity change gradually with the depth, which indicates that the entire part contains 
fresh alluvial deposits. In this part the ground water quality has played a significant role 
in resistivity data. Based on the vertical resistivity profiles it can be concluded that no 
significant contrast m the eastern side is seen in sedimentation pattern. From layered 
models a very low resistivity layer is seen which may contains salt water. This layer may 
be located near Kima at depth of 104 m, Utaria at depth 38 m with resistivity equal to 0.5 
ohm-m, Pyarepure at depth 40.5 m with resistivity 1.63 ohm-m, Bahadeona at depth 
26.28 m with resistivity 2.24 ohm-m, and also in Diwankhera located at depth 48.16 m. 
The presence of saline soil in U.P. and its relattion has been recently discussed by 
Srivastava and Singh (1998, 1999). The good fit has not been found at the locations high 
resistivity layer is present, so the thickness of layers estimated is not very accurate. The 
layered models also show that while moving from south to north along profiles I, II, EH 
and IV the resistivity values increase rapidly this may be due to the presence of fined 
sand to coarse sand to sandstone sometimes. In the western part the analysis of resistivity 
variations show sudden change in layers characteristics and their physical properties, 
which is also reflected by resistivity values. The resistivity of this part is too high 
compared to other part, which gives an indication of the presence of the layers 
characterized by different sedimentation pattern and older layers compared to other part. 
The presence of sand, compact sand or sandstone, and gravel give rise to the resistivity 
variations. The water m this part is much cleaner, this is reflected from resistivity data. 
The depth of the water table changes from this part to other part due to change in 
topography. The results show the presence of very high resistivity layer more than 1000 
ohm-m, beneath Fini^ur, Isunyian, Mahiptkhera, and Mubarkpur (Figures 4.27 to 4.30). 
This high resistivity layer is located at the depth range 40 — 80 m. From the above 
discussion, we conclude that the entire western part is uplifted area may be caused by a 
fault in the basement. Along which subsurface deformation may be taken place giving 
rise to the uplift to the basement. Such deformations may have lead to the migration of 
Ganga river towards west. 




Dcstinagar 



Figure 4.17 - App. Resi. variations (left) and Layered model estimated at Dostinagar 
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Figure 4.18 - App. Resi. variations (left) and Layered model estimated at Garah. 



Diiiakhera 



Figure 4.19 - App. Resi. variations (left) and Layered model estimated at Dullakhera. 












Figure 4.22 - App. Res. variations and Layered model estimated at Parmani 















Figure 4.24 - App. Resi. varations (left) and Layered model estimated at Kima. 



Figure 4.25 - App. Resi. variations (left) and Layered model estimated at Diwankhera. 
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Figure 4.27 - App. Resi. variations (left) and layered model 
estimated at Mubarkpur. 



Figure 4.28 - App. Resi. variations (left) and Layered model 
estimated at Muhaptikhera. 
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4.5 Integration of Remote sensing and Resistivity Data 

The resistivity data obtained from the resistivity soundings show one to one 
relation with the results obtained from the image processing (Figure 4.7). The results 
from the resistivity soundings show drastic change across a profile crossing the large 
meander cut off. Remote sensing analysis has provided a very good information about 
subsurface conditions especially when the subsurface structure cause a change to 
sedimentation pattern or any significant change to surface features like discontinuity of 
geomorpho logical features, anomalies in vegetation index, topography variations and 
change in surface layer properties. 

Thresholding techniques especially Otsu thresholding technique, which is very 
useful in geological applications and the use of different optimal thresholding values, can 
give a good idea about the study area. Thresholding has also provided information about 
lithological edges apart from the conventional edges due to roads and urban areas. 
Remote sensing data used successfully to map the fault and high land area. From the 
present study we conclude that the integration of remote sensing data and resistivity data 
is very usefiil in deducing the subsurface information and its relation with the surface 
manifestations. 
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Chapter V 


CONCLUSION 


5.1 GENERAL 

The main objective of the present study is to investigate the subsurface structure 
by using remote sensing data and to study the effect. of merging and thresholding 
techniques on geological edge extraction. To achieve these objectives the following 
studies have been carried out: 

• Study of the normalized difference vegetation index (NDVI) and classification of 
image from study area based on NDVI. 

• Map of geo morphological features from stretching FCC. 

• Various merging techniques (HPF, Price, HIS, P+XS) have been applied to improve 
the visual interpretation and to see the effect of merging on edge extraction. 

• Sobel, Laplacian, High pass filters, and Otsu thresholding technique have been 
applied to enhance the geological edges before and after merging. 

• Resistivity survey (field measurement) has been carried out from the study area. 
Resix package has been used for the analysis of resistivity soundings data. 

• The remote sensing data has been integrated with the resistivity variations obtained 
from conventional resistivity method. 

5.2 Results 

From the study carried out in present thesis the conclusion can be summarized as: 

I. Both, resistivity and remote sensing data have shown the occurring of meander 
cutoff of Ganga river. The main reasons behind occurring of channel migration 
pattern of Ganga river may be due to: 

1 . Geomorpho logical process. 

2. Occurring of cross faults in the basement. 

3. Uplift of basement rocks under the western side of meander cutoff. 
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II. Resistivity measurements have indicated a uniform change in apparent resistivity 
values on the eastern side of meander cutoff. On the other hand anonaalous 
apparent resistivity values have been found on the other side. 

III. Analysis of resistivity data has indicated that the presence of fault and high 
resistivity layer under the western side of meander cutoff. 

IV. Geological edges are very clear using Band 2 data compared to Band 3. 

V. Out of the four edge detection techniques used, Otsu thresholding technique has 
given the best edge enhancement. 

VI. NDVl contrast suggests two distinct areas which are characterized by two distinct 
subsurface settings. 

VII. The applied merging techniques have not given any improvement for geological 
edge extraction. 

VIII. Out of the four merging methods used in the present thesis the HPF and Price 
algorithm have given the best merged data. HPF is found to give much better 
results compared to Price technique. 

IX. HPF and Price algorithms have been used effectively to eliminate the cloud cover 
from merged data. 

X. Remote sensing data has proved to provide an efficient alternative to the 
conventional geological or geophysical surveys. 

5.3 Future work 

During the data collection periods, there have been a common feeling among the 
residence of that area near the fault location shows poor ground water quality and sharp 
variations in the water quality from well to well. Therefore, the study of the ground water 
quality and its relationship with surface and subsurface conditions may be taken in the 
future. 
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Appendix I 

List of tables showing the values of resistivity models for all point. 


Dostingar 


No. of 
layers 

Resistivity 
(oh nun) 

Thickness (m) 

1 

2.8 

0.935 

2 

27.18 

0.0567 

3 

31.87 

6.6 

4 

16.1 

85.12 

5 

12.1 

- 


Garaha 


So. of 
lavers 

Resistivity^ 
(ohm m) 

Thickness (m) 

1 

9.33 

0.888 

2 

14 

6.53 

3 

15.93 

28.98 

4 

24.55 

- 


Dullakhera 


No. of 
la\ ers 

Resistivity 
(ohm m) 

Thickness (m) 

1 

39.32 

0.0064 

2 

41.39 

3.97 

3 

24.58 

0.097 

4 

5.09 

4.36 

5 

44.27 

17.41 

6 

9.82 

- 


Kirna 


No. of 
layers 

Resistivity 

(ohm.m) 

Thickness (m) 

1 

83.33 

0.358 

2 

27.66 

1.74 

3 

9.82 

58.12 

4 

48.32 

43.94 

5 

0.959 

- 


Firuzpur 


No of 
layers 

Resistivity 
(ohm ni) 

Thickness (m) 

1 

297.6 

0.576 

2 

61.39 

9.08 

3 

11053 

0.591 

4 

13.72 

- 


Mubarkpur 


No. of 
layers 

Resistivity 
(ohm m) 

Thickness (m) 

1 

97.7 

0.628 

2 

55.83 

7.18 

3 

0.563 

0.2 

4 

616 

12.04 

5 

1.58 



Bhadeona 


No. of 
layers 

Resistivity 
(ohm m) 

Thickness (rn) 

1 

187.8 

1.36 

2 

132.3 

4.51 

3 

73.97 




38.13 

5 

334 

- 



































Appendix II Resistivity Observations 


Name of the place : Dostinagr 


Sr.No 

P1P2/2 

(m) 

C1C2/2 

(»”) 

Resistance 

(ohm.m) 

G 

App.Resi 

(ohm.m) 


0.5 

1 

1.4126 

2.66 

3.4 


0.5 

3 

0.254 

0.22 

7.25 

■■QHI 

0.5 

5 

0.13116 

0.08 

10.4 


1.0 

5 

0.27595 

0.16 

10.8 


1.0 

7 

0.174425 

0.0816 

13.4 

6 

1.0 

10 

0.100425 

0.04 

15.77 


1.0 

15 

0.0519 

0.0177 

18.4 



15 

0.164 


19.3 

9 



^ 0.041525 

0.0133 

19.6 

bbhobh 

5.0 

30 

1 0.07133 

0.0222 


11 

5.0 

IHESHIIi 

0.0235 

0.008 

18.45 

12 

5.0 

■HHi 

0.011525 

0.0041 

17.66 

13 

7.0 

70 

0.01683 

0.0057 

18.55 

14 

7.0 

100 

0.0091 

0.0028 

20.4 

15 

9.0 

100 

0.108 

0.0036 

18.85 

16 

11.0 

150 

0.00365 

0.00195 

11.76 

17 

15.0 

150 

0.005925 

0.00266 

13.97 


Name of the place : Garha 


Sr.No 

P1P2/2 

(m) 

C1C2/2 

(m) 

Resistance 

(ohm.m) 

G 

App.Resi 

(ohm.m) 

1 

0.5 

1 

3.885 

2.66 

9.4 

2 


3 

0.4359 



3 


5 

0.168 

0.08 

12.42 

4 

1.0 

5 

JtESSatk 

0.16 

12.86 

5 

1.0 

7 

HES^HI 

0.0816 

13.35 

6 

1.0 



0.04 

12.6 

7 

3.0 



0.12 

13.35 

8 

3.0 

15 


0.0533 

26.37 

9 

3.0 



0.0133 

14.8 


5.0 



0.0222 

14.9 

11 


■IBIHIII 

0.04 

0.0144 

17.45 

12 



0.03475 

0.0122 

17.87 






















































































Name of the place : Kirna 


Sr.No 

P1P2/2 

(m) 

C1C2/2 

(m) 

Resistance 

(ohm.m) 

G 

App.Resi 

(ohm.m) 

1 

0.5 

1 

18.6215 

2.66 

45 

2 

0.5 

3 

0.8426 

0.22 

24.06 

3 


5 

0.1756 

0.08 

16.5 

4 

■■QHH 

. 5 

0.4199 

0.16 

16.48 

5 

1.0 

7 

0.17355 

0.0816 

13.4 

6 

1.0 

10 

0.068 

0.04 

10.7 

7 

3.0 

15 

0.091 

0.0533 

WtESSaM 

8 

7.0 

30 

0-05207 

0.0311 

10.52 

9 

7.0 

50 

0.0175 

0.012 

9.82 

10 

11.0 

70 

0.015625 

0.00898 

10.9 

11 

HDEHI 

100 

0.00913 

0.0044 

13.04 

12 

hheehi 

150 

0.00545 

0.00266 

12.87 

13 

19.0 

200 

0.004 

0.0019 

13.3 

14 

19.0 

300 

0.0014 

0.00084 

mmam 

N 

fame of the place : Firozpur 

Sr.No 

P1P2/2 

(m) 

C1C2/2 

(ra) 

Resistance 

(ohm.m) 

G 

App.Resi 

(ohm.m) 

1 

0.5 

1 

84.5 

2.66 


2 

0.5 

3 

2.55 

0.22 

72.83 

3 


5 

2.25 

0.16 

88.4 

4 

1.0 

7 

0.734 


56.5 

5 

1.0 

10 

0.436 

0.04 

68.2 

6 

1.0 

15 

0.2765 


92.7 

7 

■■EHHI 

50 

0.0337 

0.0048 

44.1 

8 


70 

0.211 

0.0057 

232.6 

9 


100 

0.2291 

0.0028 

514 

10 


150 

0.0039 

0.00195 

12.56 

11 



0.0066 

0.0019 

21.1 

12 

19.0 


0.0028 

0.00084 

21 

















































































































Name of the place : Dullakbera 




17.3534 


1.506 


0.3651 


0.07653 


0.0433 


0.02115 



0.012 


0.00898 


0.0044 


41.9 


43.011 


28.7 


0.74165 

0.16 

29.12 

0.33916 

0.0816 

26 

\mmsm 

HEESflli 

25.9 

0.1733 


20.43 

0.04675 

0.0311 

22.08 


21.66 


18.9 


18.08 


12 

13.0 


70 

0.0320 

0.00266 

13 

HHEE] 


100 

0.0139 

0.0019 

14 

IfllQQ 


150 

0.0065 

0.00084 

























































































Name of the place : Diwankhera 



Name of the place : Isunyian 



































































































Name of the place : U taria 


Sr.No 

P1P2/2 

(m) 

C1C2/2 

(m) 

Resistance 

(ohm.m) 

G 

App.Resi 

(ohm.m) 

1 

0.5 

1 

11.8826 

2.66 

28.7 

2 

0.5 

3 

0.45183 

0.22 

12.9 

3 

0.5 

5 

0.12445 

0.08 

9.77 


1.0 

5 

0.2370 

0.16 

. 9.3 

5 

I.O 

7 

0.09515 

0.0816 

7.3 

6 

1.0 

10 

0.0440 

0.04 

6.9 

7 


15 

0.0696 


8.2 

8 

3.0 


0.01606 


7.58 

9 



0.01155 


6.4 

10 

7.0 

70 

0.0045 

0.0057 

4.9 

11 

11.0 

100 

0.0027 

0.0044 

3.86 

12 

15.0 

150 

o.ooi 

0.00266 

^ 2.36 

13 

19.0 

200 

0.00193 

0.0019 

6.4 


23.0 

300 

0.00074 

0.00102 

4.56 


Name of the place : Mahipatkhera 































































Name of the place : Powa 


Sr.No 

P1P2/2 

(m) 

C1C2/2 

(m) 

Resistance 

(ohm.m) 

G 

App.Resi 

(ohm.m) 

1 

IHHEHI 

1 

21.02 

2.66 

50.8 

2 


3 

0.8716 

0.22 

24.9 

3 

0.5 

5 

0.3755 

0.08 

29.5 

4 

1.0 

5 

0.7085 


27.8 

5 

1.0 

7 

0.4123 


31.75 

6 

3.0 

10 

0.667 

0.12 

35 

7 

3.0 

15 

0.25 

0.0533 

29.5 

8 

7.0 

30 

0.1716 

0.0311 

34.7 

9 

7.0 

50 

0.0518 

0.0112 

29 

10 

11.0 

70 

0.042 

0.00898 

29.4 

11 

11.0 

100 

0.018 

0.0044 

25.7 

12 

15.0 

150 

0.0109 

0.00266 

25.75 

13 

15.0 

200 

0.0057 

0.0015 

18.6 

14 

19.0 

290 

0.0025 

0.0009 

17.4 



N 

ame of the place : Pyarepur 


Sr.No 

P1P2/2 

(m) 

C1C2/2 

(m) 

Resistance 

(ohm.m) 

G 

App.Resi 

(ohm.m) 

mmmm 

■HESHH 

1 

14.174 

2.66 

34.3 

2 

jHEBHIi 

3 

1.0662 

0.22 

30.45 



, 5 

0.2688 

0.08 

21.1 


1-0 

5 

0.53125 

0.16 

20.85 


1.0 

7 

0.281 

0.0816 

21.6 

|||||[[|||||||^^ 

3.0 

10 

0.23 

0.12 

12 

HHDHH 


15 

0.11075 

0.0533 

13 



30 

0.0394 

0.0133 

18.6 ‘ 


11.0 

70 

0.09463 

0.00898 

50 

10 

11.0 

100 

0.0247 

0.0044 

313 

11 

15.0 

150 

0.01553 

0.00266 



12 

19.0 



0.0019 

19.5 

13 

19.0 



™ 0l)0084” “ 

lI.O 



































































































Name of the place : Bhadeona 


Sr.No 

P1P2/2 

(m) 

C1C2/2 

(m) 

Resistance 

(ohm.m) 

G 

App.Resi 

(ohm.m) 

1 

0.5 

1 

76.584 

2.66 

185 

2 

0.5 

3 

5.76515 

0.22 

164.6 

3 

0.5 

5 

2.06925 

0.08 

162.43 

4 

1.0 

5 

3.8665 

0.16 

151.76 

5 

1.0 

7 , 

1.53883 

0.0816 

118.5 

6 

1.0 

10 

0.859 

0.04 

135 

7 

3.0 

15 

0.4195 


49.4 

8 


30 

0.13263 


62.64 

9 

HHBHHI 

50 

0.0495 

0.0112 

27.77 



70 

0.0238 

0.00898 

16.65 

11 


100 

0.00365 

0.0044 

5.7 

12 



0.001 


2.4 

13 

19.0 

■eh 

0.00295 

0.0019 

9.75 

14 

19.0 

300 

0.00435 

0.00084 

32.5 


Name of the place : Mubarkpur 



Resistance 

(ohm.m) 

44.73 
1.9635 
0.7936 
1.7356 
0.82975 
0.3628 
0.127"~~ 
0.08485 
0.0249 
0.2745 
0.02675" " 
0.00715 
^.006 ■ 
'~0.d0165 * 

























































































APPEDIX III- Origi nal data used in this present study 
Pan, Band 1, Band 2, Band 3 













iHiiSiiSi 




,.■ • ,'y 


IRS- 1C Panchromatic Band for study area 


Band 1 (IRS- 1 CUSS HI) 



Band3(IRS-lCUSSni) 



